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Summary 


Modified barium titanate ceramics can be prepared with reduced temperature coefficients at 
room temperature. Values of coercive field and dielectric loss which permit the operation of 
transducers up to the onset of cavitation in liquids can be obtained. The mechanical and electrical 
properties of the modified ceramics change significantly with time but the drift can be minimised. 
Transducers of suitable shapes can be made by a variety of methods. 


Sommaire 


On a préparé des produits céramiques au titanate de baryum ayant de faibles coefficients de 
température 4 la température ordinaire. On peut atteindre des valeurs de champ coercitif et 
de pertes diélectriques qui permettent un fonctionnement des transducteurs jusqu’a produire 
l’apparition de la cavitation dans des liquides. Les caractéristiques mécaniques et électriques des 
produits céramiques ci-dessus varient notablement avec le temps, mais on peut réduire cette 
variation. De nombreux procédés permettent de réaliser des transducteurs de forme appropriée. 


Zusammenfassung 


Durch eine Verinderung der Zusammensetzung lassen sich Bariumtitanat-Materialien her- 
stellen, deren Temperaturkoeffizient bei Zimmertemperatur stark herabgesetzt ist. Man kann so 
giinstige Werte fiir das Koerzitiv-Feld und die dielektrischen Verluste erreichen, da die Barium- 
titanate als Fliissigkeitswandler bis zum Einsatz der Kayitation verwendbar sind. Die mechani- 
schen und elektrischen Eigenschaften weisen zwar merkliche zeitliche Schwankungen auf, die 
sich jedoch stark herabsetzen lassen. Wandler von geeigneter Form lassen sich nach verschieden- 
sten Methoden herstellen. 


1. Introduction many oxides which replace barium or titanium 
oxides without upsetting the fundamental struc- 


The characteristics ~generally-required in a 4+ of the titanate have this effect. 


piezoelectric material are a high electro-mechani- 
cal coupling coefficient and a high value of one 42000 
or other of the piezoelectric coefficients. A high ion 
coercive field, a low dielectric loss at the working 
field, small variation of properties with tempera- 


generally obtained in pure preparations of 

barium titanate and additions made for other 

purposes should therefore be kept to a minimum. 
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The general effect of additions must be consi- 
dered against the background of the properties 
of the comparatively pure ceramic material. The 
characteristic permittivity temperature curve is 
shown in Fig. 1, with transitions at 0°C and be- 
tween 100 and 120°C known as the second transi- 
tion and Curie temperature respectively. The 
other properties also change rapidly near these 
transitions. Additions may shift these transitions 
or may alter the general level of the curve. It is 
generally undesirable to lower the Curie tempera- 
ture as the material is not piezoelectric above 
this temperature. Some typical effects are listed 


in Table I. 


Table I 
Approximate effects of various 
substituents on the properties of barium titanate 
(Relative to the properties of the pure material) 


and Shape of 


A eons Effect 
re Curie point | 82" | Permuttivity | on power 
Oxide oe sition | temperature £ P 
substi- point actor 
tuted 
PbO nil 
$r0 lower nil 
Cad lower| nil | lower | broa- | lower |diminishes 
dens | peak 
MgO, MaFe |lower] nil nil | broa- | lower | diminishes 
1n0, C00, | |y sone) Pee pene 
Ni ower nil | broa- | lower | diminishes 
Na CbOs dens | peak 
Na Tas lower — broa- | lower | diminishes 
CdeCb207 dens | peak 
Y20s nil nil | lower | lower peak - 
Zr02, Sn0z | lower} lower |higher broadens nil 
ThO2 lower] lower | — broadens — 


(The words “lower peak” in the “Shape of 
permittivity temperature peak” column are 
meant to imply a lower maximum value of per- 
mittivity without any great loss in the sharpness 
of the peak: this may be attributable to the 
addition acting as a diluent, whereas when the 
peak is distinctly broadened, it must also have 
a more specific effect on the dielectric. The 40/ 
BOzg ratio is the ratio of the total number of ions 
in barium sites to the total number in titanium 
sites.) 

Lead is one of the most important substituents 
for barium. It raises the Curie point and lowers 
the second transition temperature; that is, it 
increases the relatively flat region between the 
transition points. Strontium lowers both transi- 
tion points. Calcium [1] lowers only the second— 
provided there is no excess of barium oxide. A 
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large class of small ions lower the Curie point in 
the presence of barium oxide and many of them 
have important effects on the mechanical and 
piezoelectric properties (they mostly diminish the 
piezoelectric coefficients). Yttrium [l] has a 
similar effect to calcium. Zirconia and stannic 
oxide, which replace titania, raise the second 
transition point. 

The Figs. 2 to 6 show the effects of lead and 
calcium oxides in greater detail. The results 
shown were all obtained in one laboratory with 
the same purity of raw materials and the same 
technique for obtaining a tight, internally 
strained, ceramic having a low dielectric loss at 


high fields. 


Dielectric constant %o of value at 20°C ——= 


Temperature ——= 


Fig. 2. Dielectric constant, relative to value at 20°C, 
against temperature for barium titanate modified 
with calcium and lead titanates. 

Dielectric constant at 20°C: 


1920 — Barium titanate, 

1040 — —- 5.1 weight % lead titanate, 
756 —-—-- 10.4 weight % lead titanate, 

1040 —--—-.- 6.8 weight % calcium titanate, 
538 creeeeeee 11.7 weight % lead titanate and 


a 
x44 


weight % calcium titanate. 


Fig. 2 shows the variation of the initial per- 
mittivity of polarised material. In order to show 
the variation with temperature the values have 
been plotted relative to that at 20°C. Barium 
titanate with sufficient impurities to give a tight 
ceramic structure gives a peak permittivity at 
near 0 °C and falls away steeply at lower tempera- 
tures. A ceramic containing 5 % lead titanate has 
a distinctly lower permittivity and peaks at a 
lower temperature; the variation is less than for 
the lead free material. An addition of 10% lead 
titanate further lowers the permittivity and 
gives a very low temperature coefficient at room 
temperature. 7% calcium titanate also lowers the 
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Fig. 3. Frequency constant against temperature for barium 
titanate modified with calcium and lead titanates. 
Barium titanate, 
——— 5.1 weight % lead titanate, 
—-—-— 10.4 weight % lead titanate, 
i - 6.8 weight % calcium titanate, 
seteceeeeeee 11.7 weight % lead titanate and 
7.7 weight % calcium titanate. 


&=/Y/4e 


temperature of the second transition peak [1] but 
the temperature coefficient above 10°C is quite 
high. A mixture of 12% lead titanate and 8% 
calcium titanate [2] (as given by Mason) also 
has a high coefficient above 10°C. 


Fig. 3 shows the variation of frequency con- 
stant with temperature. Barium titanate has a 
distinct minimum at 0°C which results in a 
high temperature coefficient at 20°C. The 12% 
lead, 8% calcium titanate mixture, by con- 
trast, has a minimum at some much lower tempe- 
rature and a flat maximum at 20°C giving a 
very low temperature coefficient. 7% calcium 


20 0 20 "E 40 


Fig. 4. The piezoelectric coefficient dg; against temperature 
for barium titanate modified with calcium and lead 
titanates. 

Barium titanate, 
——w— 5.1 weight % lead titanate, 
—-—-— 10.4 weight % lead titanate, 
—-—-- 6.8 weight % calcium titanate, 
seeteneenees 11.7 weight % lead titanate and 
7.7 weight % calcium titanate. 
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titanate is very little less effective than 10% 
lead while 5% lead titanate is somewhat less 
effective. 

Fig. 4 shows values of the ds; piezoelectric 
constant calculated from measurements of cou- 
pling coefficient, permittivity and frequency 
constant. The curve for the unmodified material 
shows a marked peak while the other materials 
show smaller coefficients in a similar way to their 
variation of frequency coefficient with tempera- 
ture. All the additions cause a considerable fall 
in the d coefficient. 


6-109 
Vm_}-.... 
Newton |__ 


Leal 


Piezoelectric constant g 4, —= 
e 


-20 0 20 4C 40 
Temperature —— 


Fig. 5. The piezoelectric coefficient g3; against temperature 
for barium titanate modified with calcium and lead 
titanates. 

—— Barium titanate, 
——— 5.1 weight % lead titanate, 
—:-—-— 10.4 weight % lead titanate, 
—-—-- 6.8 weight % calcium titanate, 
1.7 weight % lead titanate and 
7.7 weight % calcium titanate. 


The g3i coefficients are shown in Fig. 5. 
As these additions cause a greater fall in per- 
mittivity than in d coefficient the g values of 
the modified compositions are higher than those 
of the unmodified material at 20°C. There 
is no great difference in the temperature coeffi- 
cients. 

It can be seen from these curves that raising 
the second transition will increase the g coefficient 
as the whole temperature curve will be moved to 
higher temperatures. Zirconia and stannic oxide 
have this effect. 

Finally, in Fig.6, the electro-mechanical 
coupling coefficient is shown to be a similar 
function of temperature to the d coefficient. 
Modification with lead or calcium lowers the 
coefficient to about the same extent: the addition 
of both lead and calcium rather more so. 

Clearly no single composition is an optimum 
with respect to all the temperature characteristics 
and a suitable compromise must be accepted. 
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Fig. 6. Coupling coefficient kg; against temperature for 
barium titanate modified with calcium and lead 
titanates. 

Barium titanate, 
——— 5.1 weight % lead titanate, 
—-—-— 10.4 weight % lead titanate, 
——-+- 6.8 weight % calcium titanate, 
11.7 weight % lead titanate and 
7.7 weight % calcium titanate. 


3. Ceramic structure, coercive field and loss 


While the factors governing the temperature 
characteristics are fairly well understood, the 
control of other properties is more complex; 
many of them seem to be affected to a major 
extent by the structure of the ceramic, that is 
by the size of the crystals and by the intercrystal- 
line structure joining the crystals together and 
by the internal strains in the material. The 
coercive field may be defined as the value of the 
external field which must be applied to the 
polarised material to neutralise the effect of 
polarisation and so reduce the piezoelectric coefhi- 
cient to zero. It is not only increased by the 
substitution of lead for barium [3], which might 
be explicable on the grounds that the addition 
increases the ratio of the c to the a axis [4] at 
room temperature; but is also increased by many 
other additions which give a finely crystalline 
closely bound ceramic structure. Typical values 
are given in Table II. 


A close knit structure is present in the calcium 
lead composition devised by Mason and in this 
case it not only leads to a high coercive field but 
to a low loss at high fields. 


The general behaviour of the loss for unmodi- 
fied barium titanate is shown in Fig. 7. It rises 
sharply from the very low field value to a peak 
and then, as saturation effects become apparent, 
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Table II 
Coercive fields of modified. barium titanate ceramics 


Coercive field 
kV/cm 


Ceramic structure 


Addition 


-— Coarse crystals 2-03 

_- Glassy appearance 4.5 

5% lead titanate| Glassy appearance 9.5 

10 % lead titanate} Glassy appearance 13.0 

7% calcium Glassy appearance 7.0 
titanate 

12% lead titanate} Glassy appearance 10.0 


+ 8% calcium 
titanate 


falls off at higher fields. Even at 2 kV/cm the 
loss in this material would cause far too much 
internal heating for it to be of practical use in the 


0.4 


10 kc/s 


0.1 


0 2 4 6 
Peak a.c. field —— 


8 kV/cm 10 


Fig. 7. Variation of dielectric loss with peak alternating 
field for a typical ceramic barium titanate prepara- 
tion. 


generation of ultrasonics. In Fig. 8 another set 
of values is given for the same composition pre- 
pared under different conditions. The loss is less 
but still very high at moderate fields. The appli- 
cation of an external static field reduces the loss 
very considerably [5], but this is not a very 
acceptable way of reducing the loss in practice. 


Peak a.c. field —— 


Fig. 8. Effect of superimposed static bias on loss charae- 
teristic at 10 ke/s for a typical ceramic barium 
titanate preparation. 
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It is possible, however, by empirically adjusting 
the minor impurities and the conditions of pre- 
paration to keep the loss low up to field values 
at which cavitation occurs in most liquids. The 
combined effect of lead (10% lead titanate) and 
a carefully controlled ceramic structure is shown 


0.10 


0 4 8 12 
r.m.s. field —-> 


kV/cm 16 


Fig. 9. Variation of dielectric loss with r. m. s. alternating 
field for barium titanate containing 10% lead 
titanate. Internally strained ceramic with glassy 
fracture; 

e) 

x 


polarised, 
unpolarised. 


in Fig. 9. The power factor is kept below 5 % for 
a. c. fields up to 9 kV/cm in polarised material. 
The loss rises much more steeply with field before 
the material is polarised indicating that the 
remanent polarisation is of great importance, but 
the loss indicated by this unpolarised curve at 
2 kV/cm is still half that of the polarised material 
shown in Fig. 8 at the same field. 

Masons lead-calcium composition has a low 
loss at high fields because of the combination of 
the effects of lead and calcium on the crystal and 
ceramic structures. Only a limited amount of 
barium can be replaced by calcium; possibly 
rather less in the presence of lead; indicating that 
it is not a good fit in the barium titanate lattice 
and so will cause the strained condition which 
seems to be necessary for low loss. 

The tight, internally strained, ceramic struc- 
ture which has a low loss also has high mechanical 
strength. Coarsely crystalline structures readily 
depolarise and later shatter when repeatedly 
stressed. Fine structures sustain their properties 
under similar treatment. 


4. Ageing of barium titanate 


A property to which the approach is even more 
empirical than that of dielectric loss, is the change 
in value of the electro-mechanical properties with 
time. The permittivity, piezoelectric constants, 
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coupling factor and frequency constants of these 
materials drift by several parts in a hundred from 
their values soon after polarisation. This can be 
important in applications where the capacity or 
resonant frequency of a system is required to 
remain constant, as the drift may be sufficient to 
make the control of the temperature variation of 
the properties valueless. 

The permittivity falls even with unpolarised 
materials but Mason has argued in a recent 
paper [6] that the ageing is, nevertheless, due to 
depolarisation. 

For rough practical purposes the ageing can 
be considered to be logarithmic as 

Kk =k, + blog t. (1) 
The coefficient b may be expressed as the per- 
centage change in a parameter per decade change 
in time. 

The results on ageing are not entirely consistent 


but the following figures give the order of magni- 
tude of the effect. 


Table III 
Drift in percent per decade 


Addition Ak/k 


A®/® | Ak/K | Ad/d 


7% CaTiOs 10,5 Oo” — 2.8 — 2.5 

5% PbTiO; — 6.2 + 0.84 — 3.2 — 8.6 
k = coupling coefficient, ® = frequency constant, 
k = dielectric constant, d = piezoelectric coeffi- 
cient, 9 = density. 


As the relation between the parameters is 
k =4d©®(xo/x)” (2) 


the drift coefficients (neglecting changes in den- 
sity) should be connected by 


Ak|k = Ad/d+ A ©/O —5 Ax|k. (3) 


The drift in dielectric constant is found, as a 
rough approximation, to be proportional to the 
drift in frequency constant 


Ax/k = —AAO®/®. (4) 
A has a value of 4 for the type of internally 
strained ceramics considered here. 


Using this value and combining equation (3) 


and (4) 


neces 


Ak|k ~ Ad/d—+ Ax|x. (5) 


The general result of ageing is that the proper- 
ties of different compositions tend to converge in 
value; thus both the 7 % calcium titanate and 5 % 
lead titanate compositions have coupling coeffi- 
cients (kg;) of about 0.15 after eighteen months; 
the lead containing composition gives a 16% 
higher value when the materials are fresh. 


234 


Mason has suggested that ageing can be 
accelerated by maintaining the materials at a 
temperature of 70°C for a period of two weeks or 
so followed by a further period at an intermediate 
temperature [6]. Storing pieces at 70°C for a 
fortnight has not been found to diminish the 
subsequent ageing appreciably: but it does seem 
likely that an accelerating process could be 
devised by very slow cooling over say a month 
from a temperature near the Curie point to room 
temperature. The exact conditions would have 
to be determined experimentally to suit parti- 
cular compositions. Although means of elimina- 
ting the drift would be valuable the amount of 
ageing is not very serious for most purposes. It is 
usually more than three months after polarisation 
before a component is put into use and it will 
then be 30 months before the amount of drift 
corresponding in value to the 6 coefficients given 
in Table III, takes place. 


5. Preparation of shapes 


After devising a material with suitable pro- 
perties the final problem is to prepare it in the 
desired shape for a given application. 

Pure oxide powders do not have the plastic 
properties of clays and cannot be treated by 
conventional ceramic techniques unless suitable 
additions are made. In practice suitable additions 
can be found to make the material extrudable 
and to make a suitable slip for casting, but for 
many simple shapes a straightforward pressing 
technique gives good results. For pressing, the 
powdered material is mixed with a binder and a 
lubricant (which may be one and the same 
substance) and is pressed between punches in a 
die at pressures of 0.8 to 2.4 tons/em?. Mixtures 
of paraffin wax and stearic acid make effective 
binders and lubricants. Pressing is the most 
convenient method for bodies of uniform section 
where the pressing depth is not more than two 
or three times the die section. With careful 
technique quite long thin specimens can also be 
produced but other methods usually give better 
results in these cases. By putting in sufficient 
binder some further shaping can be carried out 
in the pressed piece: holes etc. can be drilled. 
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Where the necessity for further shaping is 
accepted hydraulic pressing can be used. In this 
case the powder is enclosed in a rubber bag, the 
bag is evacuated and is subjected to a high 
hydrostatic pressure either in air or a liquid. This 
form of pressing is more satisfactory in many 
cases than the use of punches and a die, where 
pressure is applied from only two directions. 
One difficult technical problem is the produc- 
tion of flat thin pieces of barium titanate such as 
are required for pick-up bender elements. In this 
case the process described by G. Howarr and his 
associates at the U.S. Bureau of Standards [7] 
is of value. Powdered titanate is mixed with an 
organic lacquer to make a paint and is then 
spread on a smooth surface: after the solvents 
have evaporated a thin film can be peeled from 
the base. The firing of this film is a difficult 
technical problem because of the relatively large 
area in contact with the support. The film must 
not adhere at all to its support and large thermal 
gradients must be avoided or the sheet warps. 
When removed from the furnace the sheet is 
brittle and never absolutely flat. The lack of 
precise flatness leads to difficulties if it is neces- 
sary to fix two pieces together face to face as when 
making bimorph bender elements. This difficulty 
can be avoided by painting the titanate filled 
plastic with a suspension of palladium and 
pressing the sheets together before firing [8]. 
After firing the pieces only require grinding to 
size and polarising before assembly in pick-up 


heads. 


(Received October 4th, 1955.) 
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STRUCTURAL TRANSITIONS 


IN BARIUM TITANATE-ZIRCONATE TRANSDUCER MATERIALS 


by R.C. KExu and N. J. HELLICAR 


Communication from the Staff of the Research Laboratories of The General Electric Company, Limited, Wembley, England 


Summary 


The properties of ceramic materials of the barium titanate-zirconate system having composi- 
tions with Ti : Zr ratios between 100 : 0 and 70 : 30 are described. Curves of permittivity and 
resonant frequency against temperature give a more sensitive indication of transition tempera- 
tures than X-ray powder photographs, as well as being more easily and rapidly obtained. 

The most suitable transducer material within this system is that having a Ti: Zr ratio of 
95 : 5. This material has an orthorhombic phase which is stable between 5°C and 50°C and the 
change in resonant frequency over this temperature range is only 4%. 


Sommaire 


On donne les caractéristiques de produits céramiques du type titanate-zirconate de baryum, 
ayant des rapports Ti : Zr compris entre 100 : 0 et 70 : 30. Les courbes de variation de la constante 
diélectrique et de la fréquence de résonance en fonction de la température indiquent d’une facon 
plus sensible les températures de transition que les diagrammes de poudre aux rayons X, et peu- 
vent étre relévées plus commodément et plus rapidement. 

Le matériau de ce type, qui convient le mieux pour les transducteurs a un rapport Ti : Zr 
égal 4 95 : 5. Ce matériau a une phase orthorhombique qui est stable entre 5°C et 50°C; la varia- 
tion de la fréquence de résonance dans cette gamme de températures n’est que de 4%. 


Zusammenfassung 


Die Eigenschaften yon keramischen Materialien des Barium-Titanat-Zirkonat-Systems mit 
Ti : Zr-Verhiltnissen zwischen 100: 0 und 70 : 30 werden beschrieben. Zur Ermittlung der 
Umwandlungstemperatur erhalt man durch Auftragung der Dielektrizitaétskonstante und der 
Resonanzfrequenz iiber der Temperatur ein viel empfindlicheres und auch leichter und schneller 
durchfiihrbares Verfahren als das der Rontgenuntersuchung an Bariumtitanatpulver. 

Fiir das giinstigste Material zur Verwendung als Wandler ergibt sich ein Ti : Zr-Verhialtnis 
von 95 : 5. Hierbei liegt eine orthorhombische Phase vor, die zwischen 5°C und 50°C stabil ist 


und in diesem Bereich nur eine Schwankung von 4% in der Resonanzfrequenz aufweist. 


1. Introduction 


The properties of barium titanate can be con- 
siderably modified by partial replacement of 
either or both of the cations by other metal ions 
of suitable valency, and in this way materials 
having more suitable properties for particular 
applications can be obtained. Capacitor materials, 
for example, usually require a high permittivity 
with a small temperature coefficient, but the 
material need not be ferroelectric at room tempe- 
rature. These requirements can be achieved by 
partial replacements which bring the Curie 
temperature down to room temperature and 


which produce a broad flat permittivity peak. 
The peak can be broadened if necessary by firing 
the material at a somewhat lower temperature 
than that required for a homogeneous product. 
In this way, point-to-point variations in compo- 
sition and therefore in Curie temperature are 
obtained, so that the permittivity peak will be 
broadened and flattened. 


Such piezoelectric applications as ultrasonic 
generators and detectors require a material 
having a high piezoelectric activity and a small 
temperature coefficient of resonant frequency. It 
is therefore necessary for the Curie temperature 
to be well above room temperature. It is also 
desirable for the material to be free from phase 
changes over the working range, as a change of 
direction of the polar axes occurs at the transi- 
tion temperatures. This produces partial depolari- 
sation, with a consequent reduction in piezo- 
electric activity, and rapid variations of resonant 
frequency. 

Ceramic barium titanate suffers from two 
main disadvantages: firstly it has a tetragonal 
to orthorhombic transition a little below room 
temperature, and secondly in the tetragonal 
phase the resonant frequency rises rapidly with 
temperature increasing by about 1% for every 
4°C, The investigation of the barium titanate- 
zirconate system was carried out in order to 
develop a more suitable material for transducer 
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applications. In the course of this investigation 
it also became apparent that the electrical 
measurements give a useful guide to the phase 
transitions and therefore the phase diagram can 
be produced without the necessity of taking a 
series of X-ray powder photographs at different 
temperatures on materials of various composi- 
tions. 


2. Preparation of ceramics 


Preparation follows conventional practice for 
ceramic dielectric materials. The raw materials 
are milled to ensure adequate mixing and a low 
particle size, pressed into suitable blocks and 
prefired for two hours at the temperature indi- 
cated in Table I in order to obtain the required 
chemical reactions. The blocks are crushed to 
powder, which is milled, pressed into discs of 
3.04 cm diameter und sintered by firing for two 
hours at the temperature given in Table I. An 
organic binder is added before pressing and a 
sintering agent before final firing. 


Table I 
Firing conditions and room-temperature properties 


of Ba(Ti, Zr)O, 


Firing 

temperature Room temperature results 

Ti: Zr °C 
oe Pre- | Final Radial 
fire fire coupling 
coefficient 
380) Beal! 

99 : 1 | 1250 | 1420 | 5.64 | 1350 Sok % 
98 : 2 | 1250 | 1420 | 5.64 | 1350 36.0 % 
95 : 5 | 1250 | 1420 | 5.44 | 1100 38.7 % 
92.5: 7.5 | 1250 | 1480 | 5.50 | 1350 22.0 % 
90 :10 | 1250 | 1480*] 5.54 | 1450 16.7% 
85) “215° || 1250) | 1480*|" $:55" | 2700 4 % 
80 :20 | 1250 | 1480*| 5.62 | 6750 — 
75 :25 | 1250 | 1480*| 5.24 | 5000 — 
70 :30 | 1250 | 1480*] 5.50 | 3200 — 


* Fired in oxygen atmosphere. 
£ Pp 


The compositions specified in the following 
sections indicate the proportions in which the 
raw materials are used. Although no chemical 
analysis has been carried out, it is not expected 
that any appreciable amount of any of the 
constituents will have been lost during the 
processing. 

The sintered discs are ground to a standard 
size of 2.54 em diameter and 0.254 em thick, and 
the density of each disc is determined from its 
dimensions and weight. A silvering solution is 
applied to the two major faces of the discs and 
fired on at 600 to 700°C. The material is polarised 
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for piezoelectric measurements by applying a 
field of 16 kV/cm for one hour to the dises, 


whilst they are immersed in silicone oil. 


3. Electrical measurements and results 

a) Methods of measurement 

Permittivity and loss tangent were obtained 
from bridge measurements of equivalent parallel 
capacity and resistance of the discs, the meas- 
uring voltage applied to the disc being 5 V at 
10 ke/s. Resonant frequency measurements were 
made by balancing out the parallel capacity of 
the disc and—adjusting the frequency and a 
resistance to balance a bridge at resonance. The 
electromechanical coupling coefficient was ob- 
tained from measurements of the frequencies at 
which minimum and maximum impedance occur, 
i.e. approximately the resonant and _ anti- 
resonant frequencies. All temperature runs were 
taken with the temperature rising, in order to 
avoid any temperature hysteresis effects. 


b) Room-temperature results 

Table I shows the effect of decreasing the Ti: Zr 
ratio on the electrical properties at room tempe- 
rature. The permittivity has a minimum value 
when the Ti : Zr ratio is 95 : 5 and a maximum 
value when it is 80:20. The power factor 
remains reasonably low throughout and the 
electromechanical coupling coefficient falls off 
rapidly as the replacement of titanium by zir- 
conium is increased beyond 5%. 


c) Variation of permittivity with temperature 
Permittivity curves are given in Fig. 1 for 
compositions with various Ti: Zr ratios down 
to 70 : 30 over the temperature range — 120°C 
to 120°C. The Curie temperature, at which the 
main permittivity peak occurs, decreases from 
120°C to — 60°C as the Zr : Ti ratio is changed 
from 100: 0 to 70 : 30. This fall of Curie tem- 
perature as titanium is replaced by zirconium 
accounts for the fall of coupling coefficient and 
for the rise and subsequent fall in permittivity 
shown in Table I. The permittivity peak first 


12000 


Permittivity—=> 


Ti:Zp=10:30 


~120 -80 -40 0 40 80, °C 2420 
Temperature —=— 


Fig. 1. Variation of permittivity with temperature for 
various compositions. 
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increases in height as the Ti : Zr ratio is decre- 
ased from 100: 0 to 95:5, then the peak be- 
comes lower and broader with further replace- 
ment of titanium by zirconium; this is probably 
due to small point-to-point variations in the 
material associated with inhomogeneity. A sub- 
sidiary maximum occurs in curves for materials 
having Ti: Zr above 90 : 10 at some tempera- 
ture between 0°C and 60°C. 


d) Variation of radial resonant frequency with 
temperature 


Curves of radial resonant frequency for compo- 
sitions with various Ti: Zr ratios down to 
70 : 30 are given in Fig. 2 for the temperature 
range —120°C to 80°C. These curves have been 
adjusted not to overlap excessively, by calcula- 
ting from the measured resonant frequency curve 
another curve corresponding ‘to an appropriate 
small change in diameter of the ceramic disc. 
Each curve shows two minima; no further 
minima have been observed down to —180°C. 
Taking the maximum resonant frequency bet- 
ween the two minima on the original measured 
curves for 1” discs as a criterion, the frequency 
constant is 270 ke-cm/s for BaTiO3. It is 
approximately constant at 290 ke-cm/s for 
materials having a Ti: Zr ratio between 99 : | 
and 90:10 and then rises to 350 ke-cm/s as 
the ratio is decreased to 70 : 30. 


4, Comparison of X-ray and electrical results 


X-ray powder photographs were taken on the 
materials having a Ti: Zr ratio of 95:5 over 
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Fig. 2. Variation of resonant frequency with temperature 
for various compositions. 
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the temperature range 20°C to 120°C in order 
to verify that structural changes occur at 
the same temperatures as the maxima in the 
permittivity curves (Fig. 1) and the minima in 
the resonant frequency curves (Fig. 2). These 


photographs show that there is slight depar- 


ture from cubic symmetry at 90°C, the struc- 
ture being tetragonal at lower temperatures 
down to about 50°C; at 50°C there are symp- 
toms of degeneration and below this tempera- 
ture the structure is orthorhombic down to room 
temperature. These results agree well with the 
structural conditions expected from the tem- 
perature at which the permittivity peak in Fig. 1 
and the resonant frequency minima in Fig. 2 
occur. 

Room-temperature powder photographs of 
materials with various compositions show that 
solid solution occurs over the whole range of 
compositions considered; and that a single phase 
material is produced for each composition, that 
is if the previously mentioned small point-to- 
point variations associated with inhomogeneity 
are ignored. The cell volume is increased and the 
axial ratio decreased as titanium is replaced by 
zirconium. Materials having a Ti: Zr ratio of 
97:3 and 95:5 have an orthorhombic struc- 
ture at room temperature, as would be expected 
from the curves of Fig. 2. Materials with Ti : Zr 
ratios of 92.5: 7.5 and 90:10 appear to have 
a different structure, but the diffraction lines are 
too diffuse and the departure from the cubic 
form is too small for specific identification of the 
structure. On the evidence of Fig. 2, the structure 
should be rhombohedral. Materials with lower 
Ti : Zr ratios give an essentially cubic diffraction 
pattern, so that there is no X-ray evidence that 
further replacement of titanium by zirconium 
results in the structure passing through ortho- 
rhombic and tetragonal before becoming cubic, 
as would be expected from the resonant fre- 
quency curves. However, the distortions of the 
cubic structure must be very small and the 
diffraction lines are very diffuse. 


‘ 5. Phase diagram 

The electrical results are used to build up the 
phase diagram of Fig. 3 as follows: 

(i) The Curie temperature curve is obtained 
from the permittivity peaks in Fig. 1 and marks 
the boundary between the cubic and tetragonal 
phases and therefore also between the para- 
electric and ferroelectric states. 

(ii) The other two curves are obtained from the 
resonant frequency curves of Fig. 2, The ortho- 
rhombic and rhombohedral structures are marked 
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on the diagram on the assumption that the 6. Conclusions 
minima in these curves occur at approximately (i) Solid solution occurs throughout the range 


the same temperatures as the phase changes, the of compositions investigated, in agreement with 
validity of this assumption being confirmed by results published by McQuarrie and BEHNKE[]]. 
X-ray evidence. (ii) The permittivity and resonant frequency 
curves in Figs. 1 and 2 are broadly similar to 

120 | those recently published by SMOLENSK‘ et al. [2]. | 


(iti) Electrical measurements form a reliable _ 
x and comparatively rapid means of following the 
. AA movement of the transition temperatures with — 
: etree G,. {Cubic variation in composition and are used to produce 
the phase diagram of Fig. 3. This diagram is 
a4, | confirmed by the known structural variations of 
~ barium titanate and by the X-ray evidence given 
é in this paper. 
Ss 


+ 
So 


Paar eres (iv) The most suitable transducer material 
z within this system is that having a Ti : Zr ratio 
MSR ER ae of 95:5. This material has an orthorhombic 
r phase which is stable between 5°C and 50°C and 
the change in resonant frequency over this 


temperature range is only 4%. 
(Received October 4th, 1955.) 
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THE ABSORPTION OF SOUND IN AIR AT AUDIO FREQUENCIES 
by E. J. EVANS and E. N. BAZLEY 
National Physical Laboratory, Teddington, Middlesex, England 


Summary 


The paper describes a determination of the sound absorbing properties of air at frequencies 
from 1000 c/s to 12500 c/s, for values of the relative humidity from about 5 % to 85 % at 20°C. 
The measurements were made in a large reverberation chamber, in which the surface absorption 
was very small. The absorption due to the air was evaluated by an analysis of the results of 
measurements of the empty chamber absorption. The data obtained represent a revision of 
previously published results and are in close agreement with the relaxation theory. A general 
expression is derived for the attenuation of sound in air as a function of the frequency, humidity 
and temperature. 


Sommaire 


On donne dans cet article les résultats d’une détermination de l’absorption du son dans lair 
entre 1000 et 12500 Hz, pour des taux d’humidité compris entre environ 5% et 85% a 20°C. 
Les mesures ont été faites dans une grande chambre réverbérante ot l’absorption par les parois 
était trés faible. On a déterminé l’absorption due a l’air par analyse des résultats des mesures de 
Vabsorption dans la chambre vide. Ces résultats remplacent ceux publiés antérieurement et sont 
en bon accord avec la théorie de la relaxation. On établit une expression générale de l’affaiblisse- 
ment du son dans l’air en fonction de la fréquence, de ’humidité et de la température. 


Zusammenfassung . 


Die Autoren berichten iiber die Bestimmung der schallabsorbierenden Eigenschaften von Luft 
mit relativem Feuchtigkeitsgehalt zwischen 5 % und 85 % bei 20°C im Frequenzbereich von 1000 
bis 12500 Hz. Die Messungen wurden in einem sehr groBen Hallraum mit geringer Wandabsorp- 
tion durchgefiihrt. Die gewonnenen Werte stellen eine Uberpriifung friiher veréffentlichter Ergeb- 
nisse dar und sind in guter Ubereinstimmung mit der Relaxationstheorie. Weiter wird ein allge- 
meiner Ausdruck fiir die Schallabsorption in Luft als Funktion der Frequenz, der Luftfeuchtig- 
keit und der Temperatur hergeleitet. 
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1. Introduction 


The principal features of the sound absorption 
process in air have been established by the well 
known pioneer investigations of V. O. KNuDSEN 
[1] and H. O. Kneser [2]. V.O. Knupsen by 
measurements up to 10000 c/s showed that the 
absorption of sound in air depends on the humi- 
dity and in general greatly exceeds the classical 
value for the dissipation due to shear viscosity 
and thermal conduction. Similar measurements 
in nitrogen and oxygen proved that the high 
absorption in humid air is due to the oxygen 
constituent. The results were interpreted by 
H. O. Keser on the basis of the theory of the 
molecular absorption of sound in polyatomic 
gases. He identified the anomalous air absorption 
with a relaxation process involving the vibra- 
tional energy of the oxygen molecules, and he 
accounted for the effect of humidity by taking 
the relaxation time to be a function of the water 
content. 

The measurements of V.O. KNUDSEN were 
made by the reverberation method, using two 
model chambers of different volumes to eliminate 
the effect of surface absorption. In addition 
H. Kn6rzex [3] has made measurements up to 
5000 c/s by a tube method. H. O. KnEseEr [4] 
has published a nomogram for the calculation of 
the absorption in air, based on his theory and 
KNOTZEL’s results. The values given by the 
nomogram differ appreciably from those ob- 
tained by KNuDSEN [5], and there is accordingly 
some uncertainty with regard to existing data 
and the accuracy of the theory. 

The present paper describes a determination 
carried out in the large reverberation chamber 
at the National Physical Laboratory, at frequen- 
cies from 1000 ¢/s to 12500 c/s. The work was 
undertaken as the result of the repainting of the 
chamber in 1952, when it became apparent that 
the surface absorption had been reduced to a 
remarkably low level. The sound absorbing pro- 
perties of air were derived by means of an analy- 
sis of the results of measurements of the empty 
chamber absorption as a function of humidity. 
The effects of surface, classical and molecular 
absorption in the chamber were separated by the 
procedure described below in Section 2. It was 
shown that the surface absorption in the chamber 
was small compared with the absorption due to 
the air, so that the approximate values of sur- 
face absorption furnished by this method were 
adequate for the purpose. The absorption due to 
the air was found to be in close agreement with 
the relaxation theory over a wide range of fre- 
quency and humidity. In this manner the accu- 
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racy of the theory was verified, and new and 
improved data were obtained for the absorption 
of sound in air. 


2. Theoretical considerations 
2.1. General 


The attenuation of a plane sound wave in air 

is expressed by 

I (x) = I (0) e—™2 
where I (x) = intensity at a distance x, m = inten- 
sity attenuation coefficient. 

The attenuation in air at the frequencies under 
consideration arises from the classical absorption 
due to transport phenomena, and the molecular 
absorption due to relaxation. In air the two 
effects are taken as additive and we have 

m= Me + Mp 
where m; = classical attenuation coefficient, 
Mm = molecular attenuation coefficient. 

For measurements of the decay of sound in a 
reverberation chamber, when the surface absorp- 
tion is small, the reverberation period T is given 
by the Sabine formula as modified by KNUDSEN 
to take account of air absorption, 


a: 59.3 V 1 

~ e(as+4mV) (1) 

where V = volume of chamber, a; — total sur- 

face absorption, c = velocity of sound at room 
temperature. 

The term 4mV represents the equivalent 

absorption in the room due to losses in the air. 


2.2. Classical absorption 

L. J. Stvian [6] has shown, from his own 
measurements and a survey of previous work, 
that the measured classical attenuation in air is 
of the order of 50% greater than the theoretical 
value given by the Stokes-Kirchhoff equation. 
Taking SIvIAN’s results we find for air at ¢°C 


Me = (33 + 0.2.2) f2.x 10-14 mo (2) 
where f is the frequency in c/s. 


2.3. Molecular absorption 

We adopt the theoretical treatment of H. O. 
KNeEsER [2], [7], according to which the mole- 
cular absorption in air is due to a relaxation 
process involving a lag in the adjustment of the 
vibrational energy of O2 molecules during the 
passage of a sound wave. At a frequency f we 
have to a close approximation 
t RC; Hy 


uu = Mm A = ay 2 


GG (Fe 22L) (3 


where s = m,,A = molecular attenuation per 
wave length A, k = angular relaxation frequency, 
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R=molar gas constant, C, = molar specific 
heat of air at constant volume, C; = molar 
internal specific heat due to the vibrational 
energy of the Oz constituent of air. 

The angular relaxation frequency k is approxi- 
mately the reciprocal of the relaxation time 
characterising the energy exchange between 
vibrating and non-vibrating molecules due to 
collisions. Alternatively, k is approximately the 
reaction rate constant giving the number of 
transitions per molecule per second from the 
excited to the normal state. The value of k is to 
be determined experimentally and increases with 
the humidity of the air. The water vapour acts 
as a catalyst and increases the reaction rate 
between vibrating and non-vibrating molecules. 
Previous work [1], [2], [3] indicates that the 
variation of k is small for temperature changes 
of the order of 10°C. 

On a logarithmic frequency scale, for a given 
value of k, the relation between yu and f is repre- 
sented by a characteristic bell shaped curve, 
symmetrical about a maximum at 2xf = k (see 
Fig. 3). A change in the value of k displaces the 
curve parallel to the frequency axis. The maxi- 
mum value of yw is given by 


t RC; 
Cy (Cane) (4) 


which is independent of frequency and involves 
thermal quantities only. C; is a function of 
temperature, and is obtained from spectroscopic 
data by the application of the Planck-Einstein 
relation. From the value given by HERzBERc [8] 
for the Raman displacement for Oz, the charac- 
teristic temperature for the vibration of the O2 
molecule is 2235°K. This gives for air at 20°C, 
Umax = 9.00214, and this value is adopted in the 
present work. It is in slightly better agreement 
with the experimental results than the value 
Mmax = 9.00218 given by KNEsER. 
We have further 


Umax = 


Mf 
k anf 
Es aps ©) 
where M = 2 umax/e. 


Mm = 


The value of M as a function of temperature 
is given in Fig.1. At a given frequency mn 
attains a maximum value of Mf/2 at a humidity 


such that k = 27f. 


2.4, Surface absorption 

A knowledge of the surface absorption in the 
reverberation chamber is required. We assume 
in the first instance that the chamber has ideal 
non-porous rigid surfaces, and then determine, 
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Fig. 1. Variation of M with temperature. 


as shown later, the correction to be applied for 
non-ideal surfaces and the presence of objects in 
the room. 

Absorption of sound occurs at an ideal surface 
due to viscous and thermal losses in a thin 
boundary layer of air. L. CREMER [9] has cal- 
culated the magnitude of this effect, and finds 
that the statistical absorption coefficient for 
random incidence at an ideal surface is given by 


ax =18)/f x 10-4. (6) 


2.5, Application to the reverberation chamber 

We apply the above equations for molecular, 
classical and surface absorption to the calculation 
of the absorption in the reverberation chamber 
used for the present measurements. Taking into 
account the volume and surface area of the 
chamber, the total equivalent absorption A at a 
frequency f and 20°C is, in m? units, 

ayy 137.2fx 104 | 

(7) 


+ 4.6 //f x 10-2 + 4.1 f2 x 10-8 + 6. 
The term 6 represents the allowance to be made 
for objects in the room and for the non-ideal 
surfaces. 

The value of A depends on humidity through 
the first term representing the molecular absorp- 
tion, in which k is a function of humidity. For 
the particular case k = 27f the absorption at a 
given frequency is a maximum given by 

+ 4.6) f x 10-2 + 4.1 f2 x 10-8 + 6 (8) 
which does not involve k. The correction 6 can 


accordingly be obtained from the measured value 
of Amax. Equation (7) can then be applied to the 
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evaluation of the molecular absorption at various 
humidities from the experimental determinations 


of A. 


3. Measurements 


The measurements were made in the reverbera- 
tion chamber designed for the determination of 
sound absorption coefficients at the National 
Physical Laboratory. The chamber is irregular 
in shape, with non-parallel walls and a sloping 
ceiling, and the volume is 275 m® and the surface 
area 257 m?. The walls are mainly of hard 
plaster on brickwork and the floor and ceiling of 
concrete. Steel doors 2 cm thick cover a recess 
3.1 m by 3.1 m by 18 cm deep in one of the side 
walls. Since repainting the surfaces have been 
finished in a gloss paint. 

The contents of the chamber are limited to the 
equipment used for the measurements i. e. loud- 
speakers, microphones, cables, and cords for the 

suspension of the microphones. To minimise 
sound absorption, the direct radiator loudspeaker 
units are mounted in 5 cm thick wood baffles 
on hollow concrete pillars about 1.5 m high. There 
are four loudspeakers, spaced apart on the floor 
of the chamber. 

The equivalent absorption in the room was 
determined by the reverberation method. 

Measurements were made at mean frequencies 
of 1000, 2000, 4000, 6000, 8000, 10000 and 12500 
c/s. Frequency modulated (warble) tones were 
employed, with a maximum frequency deviation 
of + 100 c/s from the mean frequency, and a 
modulation frequency of 5 c/s. 

The decay of reverberant sound in the chamber 
was registered on a calibrated high speed level 
recorder over a range of at least 60 dB. For 
measurements up to 8000 c/s four moving coil 
microphones were used, placed near the upper 
corners of the room. The method of recording 
was arranged so that a composite record was 
obtained by the superposition of the traces for 
the four microphones. The decay was substantial- 
ly linear, and the reverberation time was read 
from the average slope. At each frequency a 
number of such measurements were made and 
the mean value taken. The standard deviation 
of the individual measurements was of the order 
of 1% of the mean value. At 10000 and 12500 c/s 
a single condenser microphone, Type W.E. 640 
AA was used, placed about 2 m above the floor. 
A similar procedure was adopted, but instead of 
using different microphone positions, small 
changes were made in the mean frequency over 
a range of about 10 c/s. 

The chamber is not air conditioned, and in the 
first stage of the work results were accumulated 
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over a period of about a year as the humidity 
varied with the weather conditions. These 
measurements covered relative humidities from 
about 35% to 80%. Subsequently, the deter- 
minations were extended down to about 5% 
relative humidity by artificial drying, which was 
effected by circulating the air through drums of 
activated alumina placed in the room. The small 
amount of sound absorption due to the drying 
equipment was measured and allowed for. 

The majority of the humidity measurements 
were made with an Assman wet and dry bulb 
hygrometer, which was read in the room before 
and after the reverberation measurements. At 
very low humidities, however, this method was 
not suitable, owing to the decrease in accuracy 
of the Assman hygrometer, and because of the 
increase in humidity produced by entering the 
room. These humidities were therefore measured 
with a Meteorological Office Aircraft Frost Point 
Hygrometer, Mark II, specially designed for low 
humidity measurements, as described by Dosson, 
Brown and Cwitone [10]. This hygrometer was 
read from outside the chamber, a small quantity 
of air from the room being pumped through the 
instrument. For both hygrometers, the water 
vapour pressure in the room was obtained by the 
use of the tables of the Prussian Meteorological 
Institute [11]. The possible errors in reading the 
instruments corresponded to an uncertainty of 
the order of + 3% or less in the vapour pressure, 
the uncertainty depending on the value of the 
vapour pressure. The two instruments were 
compared at vapour pressures of about 5 mm and 
agreed to about 0.1 mm, the Assman giving the 
higher value. 


4. Results 


The results of the measurements are shown in 
Fig. 2, which gives the equivalent absorption in 
the chamber as a function of humidity at various 
frequencies at a temperature of 20°C. The curves 
drawn in the figure have been obtained by cal- 
culation as described later. 

The humidity is expressed as h, the percentage 
of water molecules in the air, i.e. the ratio of 
the number of water molecules to the total 
number of molecules present, in parts per 100, 
and given by 

10) « water jad es iad pressure i 
barometric pressure 

During the measurements the mean of the 
barometric readings was 761 mm with a standard 
deviation of 6 mm. 

The measurements were made at various room 
temperatures, with a mean value of 19°C and a 
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perimental points. Curves calculated from eq. (7) 
and (9), using the values of residual surface absorp- 
tion given in Table I. 


standard deviation of 2.1°C. The absorption 
values shown in Fig. 2 have been corrected to 
20°C by making an allowance for the variation 
of molecular absorption with temperature. The 
correction was derived with sufficient accuracy 
by making a preliminary estimate of the mole- 
cular absorption, and correcting this to 20°C by 
making allowance for the variation of M with 
temperature (Fig. 1). The measurements were 
taken over a period of more than two years, and 
the reproducibility of the results when corrected 
for temperature is seen to be satisfactory. 


5. Analysis of results 


Table I gives the maximum values of the equi- 
valent chamber absorption as calculated by eq.(8) 
on the assumption of ideal surfaces, together with 
the corresponding experimental results. 

The last column gives 6, the residual absorption 
due to non-ideal surfaces and the presence of 
objects in the chamber. The reverberation periods 
corresponding to Amax were of the order of 5 s at 


Table I 


Calculated and experimental values of Amax (m2) 


Molecular 
absorption 


Difference 
measured- 
calculated 


Classical Amax 
absorption calculated | measured 


Frequency 


c/s (max) 


absorption 
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1000 c/s and 0.5 s at 12500 c/s. The values of 6 
are approximate only, particularly at the higher _ 
frequencies, owing to the experimental error in _ 
the determination of Amax, but they are suffi- 
ciently accurate for the present purpose. The 
absorption due to the loudspeakers is known to 
be about 0.5 m2 and the remaining objects 
probably contribute a similar amount, leaving 
from about 0.5 m? to 2 m? for the residual 
absorption due to the room surfaces. The correc- 
tion for surface absorption is thus remarkably _ 
small, and it appears that the chamber approxi- — 
mates to an enclosure with ideal surfaces. The 
variation of Amax with frequency is in substantial 
agreement with theory over a wide range of 
values, and the results support the theoretical 
relations employed in the calculation. It is of 
interest to note the accuracy with which the 
values of Amax at high frequencies are predicted 
from independent data. Apart from the small 
correction to the classical absorption, the cal- 
culation involves only a knowledge of the vis- 
cosity, thermal properties and the pressure and 
density of air. 

The above values of 6 were employed to deter- 
mine the molecular absorption in the chamber by 
the application of eq. (7). From the results the 
value of w was calculated as a function of fre- 
quency at a series of values of the humidity. The 
value of k at a given humidity was then obtained 
by a curve fitting process, and a number of typi- 
cal examples are shown in Fig. 3. The bell shaped 


curve corresponds to the theoretical relation be- 
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Fig. 3. The variation of « with frequency at different 
humidities. Comparison of measured values (circles) 
with the theoretical bell shaped curve, eq. (3). 
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tween uw and f of eq. (3), and kis given by 27 fmax, 
where fmax is the frequency corresponding to 
imax: The experimental values are in very good 
agreement with theory. The values of « are more 
accurate at the higher frequencies, owing to the 
larger amount of air absorption in the chamber, 
and greater weight was given to the results above 
4000 c/s. At high humidity the results from 1000 
to 4000 c/s are seen to lie above the curves for a 
reason which is not known. As shown below, the 
effect is small when expressed in terms of mj», 
the attenuation per unit distance, but it appears 
to be quite definite. 

The variation of k with humidity is shown in 
Fig. 4 using log-log scales. The experimental 
points lie very near a straight line, and we have 
to a close approximation 


k = 1.92 hi-3 x 105 5-1. (9) 
This relation provides the information required 


for the calculation of the attenuation of sound 
in air. 


10° 


5 
005 04 


Fig. 4. The angular relaxation frequency k as a function 
of h, the per cent molar concentration of H20. 


As an overall check of the work, we apply this 
relation to the calculation by means of eq. (7) of 
the absorption in the reverberation chamber. The 
only other details required are the values of 6 
from Table I. The calculated values are given by 
the curves of Fig. 2, which are confirmed by the 
experimental points over a wide range of humi- 
dity and frequency. 

The effect previously noted at high humidities 
is again apparent in Fig. 2, more especially at 
1000 and 2000 c/s. The logarithmic scale tends 
to emphasise the discrepancy, which has a 
maximum value of the order of 0.5 m?. The 
question arises whether the additional absorp- 
tion should be attributed to the surfaces, or 
whether it is due to some secondary effect in the 
air. V.O. KNupsen [12] has found that the 
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absorption of wall surfaces is not appreciably 
affected by condensation or even by spraying 
with water. Presumably the presence of water 
would tend to decrease rather than increase the 
surface absorption, by closing up surface pores. 
It appears therefore that the increased absorption 
is probably due to a small secondary effect in the 
air. The effect does not appear to be present to 
any appreciable extent at 500 c/s. 


Relative humidity ——= 


20 30 4050 


45 10 80%. 


300 
dB/km 


200 


Intensity attenuation coefficient m ——=— 


100 


0 
0.4 06081 % 2 


h—e 


The attenuation of sound in air at 20°C as a func- 
tion of the humiditiy at different frequencies. 


Fig. 5. 


Fig. 5 gives the results obtained for the at- 
tenuation coefficient m at 20°C. The full lines 
correspond to the curves of Fig. 2 and eq. (10) 
below. The dotted portions at 1000, 2000 and 
4000 c/s are obtained from the experimental 
points of Fig. 2, on the assumption that the addi- 
tional absorption is due to the air. It is evident 
that in terms of m the uncertainty due to the 
phenomenon discussed above is small and not of 
practical significance. 

The general value of m at a temperature t°C 
and humidity h is given by 


m= Me+ My (10) 
Mf 
=f } EOS ee aia | 
m = (83+ 0.21) f2x 10+ nee 
Qnf ok 


where M is given by Fig. 1, and k by eq. (9). 


6. Comparison with previous data 


According to H. KNOrzeL the value of k lies 
between the limits given by 
k= 6.3h x 104 + 3.2 h2 x 105 


and k = 3.5 h2 x 105. 
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V. O. KNuDSEN gives k = 3.56 h2 x 10°, which 
is approximately the same as the second relation 
above. The present work gives about the same 
values for k as the above in the region of h = 0.4, 
but elsewhere the values diverge. In the case of 
KNOTZEL, however, there is better agreement if 
the experimental results are considered. His 
measurements were made at low humidities, and 
from h = 0.04 to h= 0.5 the values of k ob- 
tained are found to support the present relation. 
The discrepancy with KN6O7rzEL therefore appears 
to be attributable to the relations adopted for the 
extrapolation of the experimental data to the 
higher humidity range. 

Comparing the attenuation data, Table II 
shows that the present values differ significantly 
from the previous results. In the humidity region 
most usually encountered, from 40 % to 80 %, the 
KNESER nomogram gives consistently low values, 
and the present data lie nearer to those given by 
V. O. KnupsEn. 


Table II 
Comparison of data for air at 20°C 


Attenuation in dB/km 


Fre- 
ae Relative humidity Authority 
c/s 
10% | 20% | 40% | 60% | 80% 
1000 11 5 3 3 A 
1000 26 6 4 4 4 B 
1000 14 5 1 C 
2000 41 20 9 Te) i A 
2000 57 24 12 10 8.5 B 
2000 44 18 5 3 2 C 
4.000 110 qu 33 21 17 A 
4.000 100 84. 37 29 24 B 
4.000 lll 68 20 10 7 Cc 
6000 162 140 73 47 34 A 
6000 134. 150 69 60 51 B 
6000 140 136 49 24 16 C 
8000 196 215 126 83 60 A 
8000 155 200 103 87 73 B 
8000 157 218 84 40 28 C 
10000 220 285 191 128 94, A 
10000 185 274 141 112 98 B 
10000 180 297 120 12. 43 c 


Authorities: esent paper, 


Av Pr 
B = V. O. KnupsEn [5], 
C = H. O. Kneser [4]. 

It is of interest to consider the validity of 
eq. (10) at higher frequencies. L. J. Srvian [6] 
has published approximate data for the fre- 
quency range 15 ke/s to 500 ke/s, for air at an 
average temperature of 26.5°C and _ relative 
humidity 37%, corresponding to 1.26% H2O 
molecules. Fig. 6 gives a comparison of SIVIAN’s 
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Attenuation ——=— 


2 


‘ ie: 


Omiaae eet SaaS id 5 10° kels 
Frequency ——- 
Fig. 6. The attenuation of sound in air at 26.5°C and 37 % 
relative humidity; ---- SIVIAN, eq. (10). 


data with eq. (10). The tendency of the curves 
to approach each other above 100 ke/s is due to 
the fact that eq. (10) employs Srv1An’s value for 
the classical absorption, which is predominant in 
this region. In the anomalous absorption region 
the agreement is probably as good as can be 
expected in view of the approximate nature of 
the experimental results, and the comparison 
indicates that eq. (10) is not very wide of the 
mark at frequencies above 10 ke/s. 
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ELEKTRISCHE MODELLVERSUCHE 
ZU AKUSTISCHEN RELAXATIONSVORGANGEN 


von K. WALTHER 


III. Physikalisches Institut der Universitat Géttingen 


Zusammenfassung 


In der vorliegenden Arbeit werden die Méglichkeiten zur direkten Bestimmung von Relaxa- 
tionszeit und Relaxationsbetrag aus der Verformung eines akustischen Rechteckimpulses erértert. 
Wegen der erforderlichen groSen Bandbreite war eine experimentelle Realisierung dieser Uber- 
legungen mit den zur Verfiigung stehenden elektroakustischen Wandlern nicht méglich. Die 
akustischen Vorginge wurden im elektrischen Modellversuch nachgebildet. Das Problem der 
Trennung zweier Relaxationszeiten wurde am Modell eines elektrischen Zweipols untersucht. 


Summary 


The possibilities are considered of measuring directly relaxation time and characteristics from 
the distortion of an acoustic rectangular pulse. Because of the necessary wide bandwidth it is not 
possible to realise these considerations experimentally with the electroacoustic transducers 
available. The acoustic phenomena are translated into electrical models. The problem of the 
separation of two relaxation times is investigated with an electric two-terminal network as model. 


Sommaire 


On signale dans le présent travail différents moyens possibles pour déterminer directement 
le temps de relaxation et le taux de relaxation 4 partir de la distorsion d’une impulsion acoustique 
rectangulaire. Mais, du fait de la grande largeur de bande passante nécessaire, il n’a pas été pos- 
sible de réaliser expérimentalement ces moyens avec les transformateurs électro-acoustiques dont 
on dispose. On a représenté les phénoménes acoustiques au moyen d’un modéle électrique, et on 
a étudié sur une maquette de bipdle électrique le probléme de la séparation de deux temps de 
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relaxation. 


1. Einleitung 


Charakteristisch fiir einen thermischen Rela- 
xationsprozeB ist die endliche Einstelldauer des 
thermodynamischen Gleichgewichtes eines Teils 
der Systemfreiheitsgrade, die man ,,innere“ Frei- 
heitsgrade nennt (vgl. H. O. KneseEr [1]). Den 
sogenannten ,,auBeren‘ Freiheitsgraden eines 
Systems kann man momentan Energie zufiihren. 
Eine Energietibertragung in die inneren Frei- 
heitsgrade erfolgt iiber einen exponentiellen Aus- 
gleichsvorgang, dessen Zeitkonstante man als 
,Relaxationszeit t‘* bezeichnet. Der ,,Relaxa- 
tionsbetrag e“ stellt ein Mab fiir den Beitrag der 
inneren Freiheitsgrade zur gesamten spezifischen 
Warme des Systems dar. Im einfachsten Falle 
geniigt die Angabe je eines t- und ¢-Wertes zur 
vollstindigen Beschreibung des Relaxationsvor- 
ganges. 

Die GréBen t und ¢ kénnen aus Schallabsorp- 
tions- und -dispersionsmessungen mit sinusfér- 
miger Anregung des Systems ermittelt werden. 
Der Absorptionskoeffizient pro Wellenlange 
durchlauft als Funktion der Frequenz ein Maxi- 
mum, aus dessen Lage man die Relaxationszeit 
berechnen kann. Die Schallgeschwindigkeit steigt 


mit wachsender Frequenz um einen bestimmten 
Betrag an. Aus der Héhe der Dispersionsstufe 
erhalt man ¢. Derartige Messungen miissen in 
einem relativ groBen Frequenzbereich von etwa 
vier Oktaven ausgefiihrt werden und bedingen 
erheblichen Aufwand. 

Im folgenden soll die Frage untersucht werden, 
ob es unter Verwendung einer geeignet gewahlten 
nichtsinusférmigen Anregungsfunktion gelingt, 
die GréBen t und « eines Relaxationsvorganges 
schneller und auf direktere Weise zu gewinnen, 
als es mit der oben geschilderten Methode még- 
lich ist. Die akustischen Zustandsinderungen 
werden dabei in elektrischen Ersatzschaltungen 
modellmaBig nachgebildet. 


2. Untersuchung von Relaxationsvorgingen 
bei nichtsinusformiger Anregung 


Betrachtet man ein gasférmiges oder fliissiges 
Medium mit Relaxationseigenschaften und gibt 
an einem Ort die Schallschnelle v(t) als Anre- 
gungsfunktion vor, so laBt sich die Information 
iiber den Relaxationsvorgang aus dem zeitlichen 
Verlauf des Schalldrucks p(t) entnehmen. Fiir 
v(t) soll eine Folge von Rechteckimpulsen der 


Abb. 1. Zeitlicher Verlauf der Schallschnelle. 


Breite # und der Héhe A angesetzt werden 
(vgl. Abb. 1). Der gesamte Anregungsvorgang 
verlaufe periodisch in T. Wahlt man die Anre- 
gungsdauer @ groB gegeniiber der Relaxations- 
zeit t des Mediums und setzt T/2> 0; ¢<1 
voraus, so erhalt man mit Hilfe einer Fourier- 
Zerlegung fiir den zeitlichen Verlauf des Schall- 
druckes den Ausdruck (vgl. Abb. 2): 


See 
AD|1+ Se 


[0 
5 fir 0<t< 4%, 


E 


—A Z| 


ay zoli+te (a) 


a fir Osiss, 


P(t) = 


Dabei bedeutet Zp = 09 co den Wellenwiderstand 
des Mediums fiir die Frequenz Null. 


Abb. 2. Zeitlicher Verlauf des Schalldrucks. 


Neben dem rein ,,elastischen** Anteil im Druck- 
verlauf — einem Rechteck der Héhe Zp A — 
treten exponentielle Ausgleichsvorgange auf, wel- 
che fiir die Relaxation charakteristisch sind. Die 
GréBen t und « lassen sich, zumindest grund- 
satzlich, sehr einfach aus dem Druckverlauf p(t) 
ermitteln: t folgt aus der Zeitkonstante des expo- 
nentiellen Abklingens, ¢ ergibt sich aus dem Ver- 
haltnis der Amplitude des Ausgleichsvorganges 
zu der des elastischen Anteils. 

Will man dieses Verfahren zur Untersuchung 
von Relaxationsvorgingen experimentell reali- 
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sieren, so wird man zur Priifung der Methode 
zunachst eine Substanz mit bekanntem Relaxa- 
tionsverhalten durchmessen, bei der die GréBen 
t und ¢ méglichst hohe Werte haben sollen, um 
die Messung zu erleichtern. Bei Durchsicht der 
Literatur zeigt sich, da hierfiir nur Gase in 
Betracht kommen, und bei diesen bietet sich als 
giinstigstes Beispiel das CO g an. Bei Zimmer- 
temperatur gelten fiir Kohlensdure folgende 


Werte (vgl. [2] und [3]): 


Relaxationszeit t= 6us, 


Relaxationsbetrag «= 8%. 


Die elektroakustischen Wandler, die eine 
Impulsanregung des Gases entsprechend Abb. 1 
erzeugen oder den zeitlichen Verlauf des Schall- 
druckes wie in Abb. 2 registrieren sollen, miissen 
sehr steile Impulsflanken iibertragen kénnen. Es 
zeigt sich, daB man die Abfallzeit der Flanke 
wenigstens 10 bis 20mal kiirzer als die Relaxations- 
zeit wahlen mu, um den exponentiellen Aus- 
gleichsvorgang und damit t und ¢ noch einwand- 
frei messen zu kénnen. Bereits im giinstigsten 
Spezialfall des CO g benétigt man hierfiir eine 
elektroakustische Bandbreite von etwa 3 MHz. 
Die Anwendung eines Differenzverfahrens, bei 
dem man nur die Exponentialfunktion und nicht 
den rein elastischen Anteil im Schalldruckverlauf 
registriert, wiirde wegen der Kleinheit der Rela- 
xationsbetrage die Anforderungen an die Wand- 
ler keineswegs herabdriicken. 

Um festzustellen, welche elektroakustischen 
Bandbreiten man erreichen kann, wurden Im- 
puls-Ubertragungsversuche in Luft an elektro- 
statischen Schallwandlern mit festem Dielektri- 
kum durchgefihrt. Wirkungsweise und Aufbau 
derartiger Wandler sind in der Arbeit [4] be- 
schrieben. Als Sender und Empfanger wurden 
gleichartige Systeme verwendet, bei denen eine 
einseitig metallisierte Isolierfolie auf einer Mes- 
sing-Gegenelektrode auflag. Da eine méglichst 
hohe Eigenfrequenz erreicht werden sollte, wurde 
die Gegenelektrode hochglanzpoliert, um die 
Oberflachen-Rauhigkeit herabzudriicken und da- 
mit die Federsteife des Luftpolsters zwischen 
Folie und Metall zu erhéhen. Durch Verwendung 
einer sehr diinnen Metallpapierfolie (7 uy stark) 
wurde die schwingende Masse klein gehalten. Die 
elektrische Anregung des Sendesystems erfolgte 
durch einen StoB-Oszillator mit einer Schaltung ~ 
analog wie in [5], S. 207. Die Abb. 3 zeigt den 
elektrischen Sendeimpuls. Auf der Empfangs- 
seite wurde eine méglichst hohe elektrische Vor- 
spannung an das System gelegt, um die Feder- 
steife des Luftpolsters und damit die Eigen- 
frequenz zu vergréBern. Der Empfangsimpuls 
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Abb. 3. Elektrischer Abb. 4. Akustischer Impuls am 


Impuls am Sendesy- Empfangssystem; Ausschwing- 
stem; Impulsbreite 0.8 periode 2,4 wus, Impulshéhe 
us, Impulshéhe 110 V. 80 mVgg. 


bei einer Vorspannung von 280 V ist in Abb. 4 
dargestellt. Man erkennt, daB diese elektrostati- 
schen Wandler nicht in der Lage sind, die steilen 
Flanken des elektrischen Anregungsimpulses zu 
iibertragen, sondern nur einen stark gedimpften 
Ausschwingvorgang in ihrer Eigenfrequenz von 
etwa 420kHz zeigen. Diese Ubertragungsver- 
suche lieferten also nicht die erforderliche elektro- 
akustische Bandbreite von einigen MHz. 

Die Verwendung eines Wandlers, bei dem ein 
piezoelektrischer Schwinger (BaTiO3) einseitig 
an ein Material vom gleichen Wellenwiderstand 
(Messing) angekoppelt wurde, lieB es aussichts- 
reich erscheinen, daB man in gasférmigen Medien 
Ubertragungsbandbreiten von einigen MHz er- 
reichen kénnte. Doch zeigte eine quantitative 
Abschatzung, da die zu erwartenden elektri- 
schen Empfangsspannungen bei Impuls-Uber- 
tragungsversuchen sehr gering sind und bei der 
groBen Bandbreite bereits im thermischen Rau- 


schen liegen. 
Zusammenfassend laBt sich sagen: Bei Unter- 
suchung von Relaxationsvorgingen mit sinus- 
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3. Nachbildung akustischer Relaxationsvorginge 
im elektrischen Modellversuch 


Bedient man sich einer Analogie zwischen 
akustischen und elektrischen GréBen, so lassen 
sich die akustischen Zustandsanderungen eines 
relaxierenden Systems am Modell eines Relaxa- 
tionszweipols elektrisch nachbilden. Eine még- 
liche Ersatzschaltung ist in Abb. 5 
stellt. (U=Spannung, [=Strom). Die Konden- 
satoren C’ und C bilden die Energiespeicher 
fiir die ,,auBeren“ und ,,inneren“ Freiheitsgrade 
des Systems nach. Man erkennt unmittelbar, 
daf eine Energiezufuhr an die auBeren Freiheits- 


darge- 


grade momentan erfolgen kann, wihrend die 
Energieiibertragung auf die inneren Freiheits- 
grade mit einer endlichen Einstellzeit vor sich 
geht. Ein derartiges Modell bildet, akustisch ge- 
sehen, die Zustandsanderungen in einem Volum- 
element nach, dessen Ausdehnung klein zur 
Schallwellenlange ist. Die Zustandsgleichung fiir 
ein kompressibles Medium mit einer Relaxa- 


tionszeit lautet: (vgl. [13], S. 366) 
' i) Hare , do 
(a ne =a(e+e z}. (2) 


Dabei bedeuten: 
p = Druckanderung, 9 = Dichteanderung, 
Cy = Schallgeschwindigkeit fiir die Frequenz 
Null, t’ = ,,Relaxationszeit**, tT” = ,,Retar- 
dationszeit“‘; t’ und +” sind fiir kleine Re- 
laxationsbetrage annahernd gleich. Die Differen- 
tialgleichung des elelektrischen Zweipols nach 
Abb. 5 hat eine analoge Form wie Gl. (2), wenn 
man p durch U und oe durch Q = {Ide ersetzt. 


Abb. 5. 
Relaxationszweipol mit 
einer Zeitkonstante. 


formiger Anregung mu man Relaxationszeit und 
Relaxationsbetrag aus einer gréBeren Anzahl von 
Einzelmessungen ermitteln, wobei die Bandbreite 
bei jeder Messung grundsatzlich sehr klein sein 
kann. Bei Anwendung von nichtsinusférmigen 
Anregungsfunktionen kénnte man die GréBen tr 
und ¢ direkt aus einer einzigen Messung heraus 
erhalten. Wegen der erforderlichen groBen Band- 
breite war eine experimentelle Realisierung dieses 
Verfahrens mit den zur Verfiigung stehenden 
elektroakustischen Wandlern nicht méglich. 


Abb. 6. 
Glied des Kettenleiters 
ohne Relaxationsverluste. 


aber si 
=o) 
Abb. 7. 


Glied des Kettenleiters 
mit Relaxationsverlusten. 


Will man den Vorgang der eindimensionalen 
Schallausbreitung im ralaxierenden Medium 
durch einen elektrischen Modellversuch nach- 
bilden, so muf man von Zweipolen mit kon- 
zentrierten Schaltelementen zu Vierpolen mit 
,,verteilten® elektrischen Parametern iibergehen, 
d.h. man muf eine homogene elektrische Lei- 
tung mit Relaxationsverlusten betrachten. Die 
Wellenausbreitung auf einer relaxierenden Lei- 
tung wird durch folgende Differentialgleichung 
beschrieben: 
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0 (ew pune ey 50 uw 
. = | — ¢, : 3 
"ot ( Openiae = ( ae Oe (3) 


Dabei bedeuten: 
Cy = Phasengeschwindigkeit fiir die Fre- 
quenz Null, c.. = Phasengeschwindigkeit 


fiir die Frequenz Unendlich, t = Relaxa- 


Coo 


2 
—€¢ 

0 e 

— = Relaxationsbe- 


tionszeit, ¢ = 5 
trag. “co 
Die Funktion y= y(x,t) kann Strom oder 
Spannung am Ort x der Leitung zur Zeit t dar- 
stellen. Die entsprechende akustische Wellen- 


gleichung erhalt man, wenn man mit Hilfe der 


Druck p —— Spannung U, 
Schnelle v —- Strom I 


Analogie die elektrischen durch die akustischen 
GréBen ersetzt. 

Eine derartige elektrische Leitung wurde durch 
einen Kettenleiter aus konzentrierten Schalt- 
elementen angendhert, weil die Realisierung eines 
homogenen Kabels mit relaxierendem Dielek- 
trikum im Bereiche tiefer Frequenzen schwierig 
ist. Diese Annaherung ist nur im Durchlaf- 
bereich des Vierpols fiir Frequenzen weit unter- 
halb der Grenzfrequenz statthaft. Zur Verfiigung 
standen zwei Laufzeitketten, die aus je 40 
T-Gliedern in der Schaltung nach Abb. 6 be- 


standen. In eine Laufzeitkette wurden kiinst- 
lich Relaxationsverluste durch Einbau von RC- 
Gliedern parallel zu C’ eingefiigt (vgl. Abb. 7). 
Die Abb. 8 und 9 zeigen die Ausgangsspannungen 
bei Ubertragung eines elektrischen Rechteck- 
impulses durch beide Laufzeitketten. Der Ketten- 
Jeiter ohne Relaxationsverluste iibertragt den 
Impuls praktisch unverzerrt. Die kleinen Ein- 
schwingvorginge hinter den vertikalen Flanken 
sind durch die TiefpaScharakteristik des Vier- 
pols mit einer Grenzfrequenz von etwa 20 kHz 
bedingt. Beim Kettenleiter mit Relaxationsver- 
lusten zeigt der Spannungsverlauf im Anschlub 
an die Impulsflanken exponentielle Ausgleichs- 
vorgange. Man kann die GréBen +t und « des 
Relaxationsvorganges unmittelbar aus dem Oszil- 
logramm entnehmen. 


Abb. 8. Impuls am Ausgang 
des Kettenleiters ohne Rela- 
xationsverluste. 


Abb. 9. Impuls am Ausgang 
des Kettenleiters mit Rela- 
xationsverlusten. 
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4, Elektrische Modellversuche zur Trennung 


zweier Relaxationszeiten 


Die bisherigen Betrachtungen bezogen sich auf 
Relaxationsvorgange, die durch Angabe einer 
Relaxationszeit und eines Relaxationsbetrages 
vollstandig gekennzeichnet werden. Im allge- 
meinen muB man jedoch eine Uberlagerung sehr 
vieler Relaxationskomponenten zur Darstellung 
der MefBergebnisse ansetzen. Die rechnerische 
Analyse eines Relaxationsspektrums aus Mef- 
werten bei sinusférmiger Anregung des Systems 
behandeln z. B. die Arbeiten [6], [7], [8] und [9]. 
Der allgemeine Fall ist fir eine Nachbildung im 
elektrischen Modellversuch viel zu kompliziert. 
Es soll hier nur der einfachste Spezialfall eines 
diskreten Spektrums aus zwei Relaxationszeiten 
betrachtet werden, wie er z. B. bei der Kohlen- 
sdure vorliegt. Hier kénnen zwei verschiedene 
Normalschwingungen des CO2-Molekiils als innere 
Freiheitsgrade angeregt werden (vgl. [10] und 
(uu). 

Fir den elektrischen Modellversuch wurde 
eine spezielle nichtsinusférmige Anregungsfunk- 
tion verwendet, welche bis zu einem gewissen 
Grade eine Trennung zweier Relaxationszeiten 
erméglicht. Der betrachtete Zweipol mit zwei 
Zeitkonstanten ist in Abb. 10 dargestellt. Als 


Tit) 


Abb. 10. Relaxationszweipol vi) 
mit zwei Zeitkonstanten. 


Anregungsfunktion wurde die Spannung U/(t) 
vorgegeben. Die Information iiber den Relaxa- 
tionsvorgang ist dann im Strom I(t) enthalten. 
U(t) wurde als ,,stiickweise lineare Zeitfunktion™ 
angesetzt, die sich mit der Periode 27 wiederholt. 
Thr Verlauf ist in Abb. 11 skizziert. Im Zeitinter- 


Abb. 11. Zeitlicher Verlauf der Anregungsfunktion. 


vall t = 0---t erfolgt ein linearer Anstieg um den 
Betrag A, dann ein Sprung der Héhe 2Up, daran 
schlieBt sich im Bereich t = t--:27 ein linearer 
Abfall um den gleichen Betrag A an. Analytisch 
ist U(t) durch folgenden Ausdruck gegeben: 
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(42 +3) anes 
U(t) = a: (4) 
—(4 —* +B) fare en Shes 2T 
A 
eA O 5 - 


Der Strom J(t) im eingeschwungenen Zustand 
wurde nach einem Verfahren von D. L. WarpE- 
Licu [12] berechnet. Man erhalt als Ergebnis: 


1. Fir das Zeitintervall 0 <t <tT: 


mere OY (e%—e?%) 
i Wee é [ite tet zy) 
2 @e22, 
tea —eray (aed) em 


2. Fur das Zeitintervall +t <t < 2r: 


A Cc ye gue Dy 
I(t) =— [i+ éeytegt 2e1 even 
>< (= G1 =1) e—(t—aDitn + (6) 
re) (Oo 5, EM 
+ 2 €9 (em) (a xea—lje ‘ 


Die Funktion verlauft periodisch in 27. Dabei 
bedeuten: 


tT = RC, te = ReCo die Relaxationszeiten, 

€) = C1/C’, e2 = Ce/C’ die Relaxationsbetrage; 

xy = t/t1, %X%2—=T/t2 sind dimensionslose 
Zahlen. 

Man erkennt aus Gl. (5) und (6), daB der Strom 
durch den Zweipol neben einer Rechteckkompo- 
nente noch zwei exponentielle Ausgleichsvor- 
ginge mit den Zeitkonstanten tT; und T2 enthalt, 
deren Amplituden, bezogen auf den ,,rein ela- 
stischen“ Anteil, gréBenordnungsmabig durch 
die Relaxationsbetrage 2; und eg gegeben sind. 

Das Blockschema des elektrischen Modellver- 
suchs zur Trennung zweier Relaxationszeiten ist 


I 


: 


C’ 
1 |G 
Zweipol 1 |_ —7— 


Blockschema des Modellversuches. 
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in Abb. 12 skizziert. Es wurde ein Differenzver- 
fahren angewendet. Der Impulsgenerator zur 
Erzeugung der Spannung U(t) hat in der End- 
stufe zwei Ausgange, an denen die erzeugte 
Impulsform mit entgegengesetzter Polaritat ab- 
genommen werden kann. Die beiden Zweipole I 
und II sind gegeneinander geschaltet. Die Diffe- 
renz der Stréme durch I und II wird an einem 
sehr kleinen Widerstand Rs in einen proportio- 
nalen Spannungsabfall umgesetzt, der auf einem 
Kathodenstrahl-Oszillographen angezeigt wird. 
In der akustischen Analogie kann z. B. der Zwei- 
pol I eine mit COg gefiillte Druckkammer und 
der Zweipol II eine entsprechende mit Luft ge- 
fiillte Kammer darstellen. Durch dieses Differenz- 
verfahren kann man den Relaxationsanteil vom 
,rein elastischen’* Anteil trennen und auf diese 
Weise die exponentiellen Ausgleichsvorgiange, 
welche die Information iiber die Relaxations- 
zeiten enthalten, gut beobachten. 

Als Vorstufe sei zunadchst der Fall diskutiert, 
daf der Zweipol I nur eine Relaxationszeit 71 
enthalt. Gibt man im Anregungsvorgang U/(t) die 
Periode 2t sowie den Anstieg A fest vor und 
variiert die Sprunghéhe 2U9 mit Hilfe des 
Dampfungsgliedes (vgl. Abb. 12), so erkennt man 
aus Gl. (5), daB bei Kinstellung des Wertes 


Uo/A = 11/7 (7) 


die Exponentialfunktion mit der Zeitkonstante 11 
verschwindet. Der Strom durch den Zweipol ver- 
Jauft dann rein rechteckférmig. Der Kiirze halber 
soll fiir diesen Zustand die Ausdrucksweise ,,Die 
Zeitkonstante t, ist kompensiert worden“ ver- 
wendet werden. Ist t bekannt und mit man 
Up/A aus dem Oszillogramm von U(t) nach 
Abb. 11, so bietet Gl. (7) die Méglichkeit, eine 
unbekannte Relaxationszeit tT; zu berechnen: 
mesa Ly Vee (8) 
Aus Gl. (5) geht ferner hervor, daB fiir Up/A < t,/t 
baw. Uo/A > 1%/t die zur Zeit t = 0 einsetzen- 
den exponentiellen Ausgleichsvorgainge negati- 
ves bzw. positives Vorzeichen haben. Diese Zu- 
stande seien kurz mit ,,Unterkompensation~ 
bzw. ,, Uberkompensation“* der Zeitkonstante Tt, 
bezeichnet. Die Spannungs- und Strom-Oszillo- 
gramme der Abb. 13 und 14 sollen diese Verhalt- 
nisse verdeutlichen. 
Geht man anschlieBend zu dem Fall tiber, daB 
der Zweipol I zwei Relaxationszeiten 1; und t2 
enthalt, so erkennt man aus Gl. (5), dafs man 


durch die Vorgabe 
(Uo/A)i = t/t (9) 
die Zeitkonstante tT; und durch 


(Uo/A)2 = T2/T (10) 
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Abb. 13. Spannung U(t) Abb. 14, Strom I(t) 
Zweipol mit einer Relaxationszeit ; 
(a) Unterkompensation, (b) Kompensation, 
(c) Uberkompensation. 
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schiedliche Zeitkonstanten im Zweipol I ist der 
Verlauf des Stromes I(t) in Abb. 15 dargestellt, 
wobei in der Reihenfolge dieser Bilder die Sprung- 
héhe Up/A in der Anregungsspannung U(t) zu- 
nimmt. Die ,,kurze“‘ Zeitkonstante des Zweipols 
sei mit tT. bezeichnet. Sie wird bereits bei einer 
kleineren Sprunghéhe kompensiert (Abb. 15b) 
als die ,,lange“‘ Zeitkonstante 1; (Abb. 15e). 

Die Trennung der beiden Komponenten des 
diskreten Relaxationsspektrums ist in Abb. 15 
sehr gut méglich, weil hier die Verhaltnisse fiir 
den Modellversuch besonders giinstig gewahlt 
wurden. Eine  genauere graphische Diskussion 
der Impulsformen nach Gl. (5) und (6) zeigt, 
daB eine getrennte Messung zweier Relaxations- 
komponenten mit gleichen Relaxationsbetragen 
nicht mehr mit Sicherheit méglich ist, wenn die 
zugehérigen Zeitkonstanten sich um wesentlich 
weniger als etwa einen Faktor 5 voneinander 
unterscheiden. In dieser Hinsicht vermag das 
beschriebene Analyseverfahren mit nichtsinus- 


Abb. 15. Strom durch einen Zweipol mit zwei Zeitkonstanten; 


(a) 71, T2 unterkompensiert, 

(b) t1 unterkompensiert, t2 kompensiert, 

(c) tT, unterkompensiert, Tz tiberkompensiert, 
(d) t, unterkompensiert, Tz tiberkompensiert, 
(e) t1 kompensiert, t2 tiberkompensiert, 

(f) 71, T2 tiberkompensiert. 


kann. 
».Xompensation~ bedeutet in diesem Falle, daB 
jeweils nur ein exponentieller Ausgleichsvorgang 
auftritt. Man kénnte diese Zustainde etwa durch 
Aufzeichnung iiber einen logarithmischen Ver- 
starker feststellen und die zugehérigen Zeit- 
konstanten aus Gl. (9) und (10) berechnen. Aller- 
dings mi®te dann der Rechteckanteil in der 
Differenz der Stréme durch die Zweipole I und II 
véllig kompensiert werden, was in Abb. 15 noch 
nicht geschehen ist. Bei einem Spektrum mit 
mehr als zwei Relaxationszeiten versagt dieses 


die Relaxationszeit t2 ,,kompensieren“ 


Verfahren grundsatzlich. Fiir zwei sehr unter- 


formiger Anregungsfunktion nicht mehr zu lei- 
sten als die Untersuchungsmethoden mit sinus- 
formiger Anregung. Haben die Relaxations- 
betrage ¢1, ¢2 beider Komponenten sehr unter- 
schiedliche GréBen, so wird die Analyse des Vor- 
ganges mit dem kleineren Betrag noch wesentlich 
problematischer als im oben diskutierten Falle. 

Herrn Prof. Dr. E. Meyer danke ich fir die 
Anregung zu dieser Arbeit und fiir sein stets fér- 
derndes Interesse. Ebenso sage ich den Herren 
Dr. K. Tamm Dr. G. Kurtze herzlichen 
Dank fiir viele wertvolle Diskussionen. 


und 


(Eingegangen am 10. September 1955.) 
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LE THERMOPHONE, 
PISTONPHONE AERODYNAMIQUE 


par P. Rrery 


Ancien éléve de I’Ecole Polytechnique, Centre National d’Etudes des Télécommunications, Paris 


Sommaire 


On élabore ici une nouvelle théorie du thermophone basée sur l’application de principes 
généraux de la thermodynamique et faisant appel 4 V’interprétation des phénoménes physiques. 
On calcule la pression acoustique produite 4 l’intérieur de la cavité du thermophone. On retrouve 
une formule établie précédemment par l’auteur a la suite de considérations théoriques rigoureuses. 
On espére que cette nouvelle formule donnera des résultats en bon accord avec ceux déduits des 
autres méthodes utilisées actuellement pour l’étalonnage des microphones. 


Zusammenfassung 


Ausgehend von allgemeinen thermodynamischen Prinzipien gibt der Autor eine neue Theorie 
des Thermophons. Aus der Berechnung des im Innern des Thermophon-Hohlraumes erzeugten 
Schalldrucks ergibt sich eine Formel, die vom Autor schon vor kurzem auf Grund sorgfaltiger 
theoretischer Uberlegungen aufgestellt worden war. Man hofft, daB man mit Hilfe dieser neuen 
Formel gute Ubereinstimmung mit anderen Methoden zur Mikrophoneichung erzielen kann. 


Summary 

A new theory of the thermophone is worked out based on the application of general thermo- 
dynamic principles and applied to the interpretation of the physical phenomena. The acoustic 
pressure inside the thermophone is calculated anew and found to agree with an earlier formula 
of the author’s derived by rigorous theory. The new formula should give results in accordance 


with those of the methods commonly used for calibrating microphones. 


Avant-propos 


Nous avons exposé, dans un récent numéro des 
Annales des Télécommunications [1], une théorie 
nouvelle du Thermophone, qui, tout en s’inspi- 
rant des théories classiques, présente le calcul de 
la pression produite dans une cavité par un ther- 
mophone a feuilles d’or, de maniére plus rigou- 
reuse et plus compléte. Nous voudrions ici re- 
prendre le probléme sous un jour différent, en 
insistant plus particuliérement sur laspect 
physique des phénoménes, car le thermophone 
est, non seulement un appareil qui a rendu de 
grands services en acoustique physiologique ou 


pour l’étalonnage des microphones de référence, 
mais aussi un instrument de physique parti- 
culiérement intéressant, puisqu’il permet d’émet- 
tre un son effectif sans l’intermédiaire d’aucun 
organe mécanique vibrant. 


Introduction 


Lorsqu’on examine les théories successives qui 
ont été élaborées sur le thermophone durant ces 
trente ou quarante derniéres années [2]---[6], on 
est frappé par deux choses. Tout d’abord il 
semble, qu’a part peut étre ARNOLD et CRANDALL 
[2], les différents auteurs ne se soient pas pré- 
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occupés de l’aspect physique des phénoménes. 
Ensuite, et c’est peut étre la conséquence de la 
premiére remarque, on constate que toutes les 
théories sont calquées sur le méme processus de 
calcul, chaque auteur s’étant contenté de modi- 
fier quelques hypothéses secondaires. 

Or ce processus souléve de grandes difficultés. 
Il nécessite en effet la détermination de létat 
d’un gaz en tous les points d’une cavité. Dans le 
cas présent, comme il s’agit d’une cavité irré- 
guliére, on est amené a faire des approximations 
souvent trés grossiéres et peu justifiées, ce qui a 
entrainé de trés vives critiques de la théorie. De 
plus les auteurs paraissent avoir attaché une trop 
grande importance aux feuilles d’or qui jouent 
dans leurs calculs un réle presque aussi consi- 
dérable que la cavité, alors que nous verrons 
qu’il est beaucoup plus logique, dans l’interpréta- 
tion des phénoménes, de considérer les feuilles 
comme un élément extérieur au systéme constitué 
par la cavité, élément perturbateur qui absorbe 
une partie de l’énergie thermique apportée au 
gaz. On pourrait parfaitement concevoir, par 
exemple, un thermophone a gaz ou l’énergie 
calorifique serait apportée par une flamme dont 
on ferait varier périodiquement la longueur.! 
Nous traiterons donc d’abord le probléme du 
thermophone en faisant abstraction de la source 
physique de chaleur; et nous serons ramenés a 
Vétude d’un gaz enfermé dans une cavité, en 
équilibre mécanique, mais non en équilibre 
thermique. Le comportement de ce gaz s’explique 
alors par une succession de transformations éner- 
gétiques simples, ce qui permet, par une appli- 
cation immédiate des principes généraux de la 
thermodynamique, de calculer la pression acous- 
tique engendrée dans la cavité. 

Ce nouveau processus offre le double avantage 
d’étre rigoureux du point de vue théorique, et de 
donner une vue d’ensemble sur l’aspect physique 
des phénoménes. On peut alors faire intervenir 
Pinfluence des feuilles d’or, ainsi que celle de 
certains phénoménes secondaires, qui se traduit 
sous la forme de termes correctifs apportés dans 
la formule donnant la valeur de la pression acous- 
tique. 

Nous serons amenés, pour mieux dégager les 
phénoménes ‘essentiels, 4 admettre des hypo- 
théses simplificatrices au cours de cet exposé. Les 
calculs seront done approchés, mais ce qui im- 
porte ici c’est le raisonnement physique, toujours 
un peu schématique, et qui parle beaucoup plus 
a Pesprit qu’une formule rigoureuse mais compli- 
quée. La théorie compléte que nous avons éla- 


1 Ce thermophone & gaz se rapprocherait de la flamme 
chantante classique. 
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borée d’autre part nous conduit d’ailleurs 4 des 
résultats suffisamment semblables pour que les 
simplifications admises ici puissent étre consi- 
dérées comme physiquement valables. 


1. Théorie élémentaire du thermophone 


Nous considérerons seulement ici, pour bien 
délimiter le probléme, le cas du thermophone a 
feuilles, d’or en atmosphére d’hydrogéne. Il est 
d’ailleurs 4 noter que l’étude du comportement 
du gaz reste identique dans le cas du thermo- 
phone 4 fil, et qu’elle ‘subsiste méme quelle que 
soit Vorigine de l’énergie calorifique apportée au 
systéme, qu’elle soit électrique, chimique, ou 
autre. 

Le thermophone considéré est constitué par 
une cavité étanche et rigide en laiton, ayant la 
forme d’un cylindre aplati, et remplie d’hydro- 
gene. Deux feuilles d’or de 1/10 de micron 
d’épaisseur sont tendues sensiblement au milieu 
de cette cavité entre des plots également en laiton. 
Le schéma d’ensemble est donné par la Fig. 1. 


Feuilles 


ZH 


a7 
PE NIZ 


Plot Feuilles 


Fig. 1. Vue schématique du thermophone. 


Le fonctionnement du thermophone est alors 
le suivant: Les feuilles d’or sont parcourues par 
un courant périodique de fréquence f. Elles 
s’échauffent et se refroidissent périodiquement 
par effet Joule, et font varier par conduction la 
température du gaz environnant. Ces variations 
thermiques entrainent au sein du gaz des varia- 
tions de volume, et par conséquent des variations 
de pression, c’est-a-dire une pression acoustique. 
Nous avons done affaire 4 un volume gazeux 
enfermé dans une cavité rigide, et il s’agit de 
déterminer la pression acoustique produite par 
une variation de l’état thermodynamique du gaz, 
celui-ci échangeant périodiquement avec l’exté- 
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rieur, par l’intermédiaire des feuilles d’or, une 
certaine quantité d’énergie sous forme calorifique. 

On voit tout de suite, ce qui sera précisé dans 
la suite, une certaine analogie entre le thermo- 
phone et le pistonphone: dans les deux cas nous 
retrouvons le volume gazeux, enfermé dans une 
cavité, et dont l’état thermodynamique varie. La 
seule différence est que, dans le cas du piston- 
phone, l’échange d’énergie avec l’extérieur s’effec- 
tue au moyen du travail d’un piston, alors que, 
dans le cas du thermophone, il s’effectue sous 
forme calorifique. Nous allons voir que les calculs 
sont analogues dans les deux cas et peuvent se 
traiter parallélement. 

Considérons le systéme formé par le volume 
gazeux V enfermé dans la cavité. Nous suppose- 
rons que le gaz se comporte comme un gaz parfait 
et nous supposerons aussi que la variation de]’état 
thermodynamique du gaz s’effectue suffisamment 
lentement pour que le gaz se trouve toujours trés 
prés dune position d’équilibre mécanique. Cela 
signifie que la pression instantanée du gaz est prati- 
quement la méme a chaque instant en tous les 
points de la cavité. Cette hypothése se trouve 
vérifiée lorsque la fréquence de la variation est trés 
inférieure a toute fréquence de résonance de la 
cavité. Dans ces conditions, la transformation de 
Vétat du gaz est mécaniquement réversible. Elle 
est parfaitement déterminée si on connait l’état 
initial et état final. Cet état dépend de deux 
variables indépendantes, la pression et le volume, 
par exemple. 


Y YY 
Fig. 2. Schéma du pistonphone. 


a) Cas du pistonphone (Fig. 2) 

Nous allons considérer deux états voisins du 
gaz: Le premier est défini par une pression p, un 
volume V, le deuxiéme par une pression p+ dp, 
un volume V—dv. Le systéme a recu du milieu 
extérieur le travail dr—pdv mais il n’a pas 
échangé de chaleur avec lui. Or, la quantité de 
chaleur dQ échangée au cours d’une transforma- 
tion réversible s’exprime, en fonction des deux 
variables choisies, sous la forme: 


dQ = Adv + udp. (1) 
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Les facteurs / et yu se calculent?, connaissant les 
capacités calorifiques du gaz a pression constante 
Cp, et a volume constant cy: 


{SEE SON Rae) 
Ais \on/Ngniee Re 

Will 5000s wit oer a 
ee (inion on 


puisque p, V et T sont liés par la relation des gaz 
parfaits: 


paVi=sheT 


Si on exprime que dQ est nul, on trouve la 
relation entre dv et dp: 


1 Vet A adi Ve Pip de, 
dp = F dv = V 
ou dp = y dE/V, (2) 


avec y = Cy/c, et dE énergie apportée par le 
milieu extérieur. 


b) Cas du thermophone 

Nous allons de méme considérer deux états 
voisins du gaz. Dans le cas présent, le volume n’a 
pas varié, le travail fourni par l’extérieur est 
done nul. Mais le milieu extérieur a fourni une 
quantité de chaleur dQ, qui s’exprime, comme 
précédemment: 


dQ = Adv + yw dp. 
Mais ici, dv est nul. On obtient donc: 


nd Oy, ed Ores Re ak, 
er i Dad AlN be eal i fi 


dp 


ou encore, en utilisant la relation de Mayer 


Cp — Cy = (v—1) m= RJ: 
dp = (y—1) dk/V. (3) 


2 Le calcul de A et wz s’obtient en s’appuyant sur les pro- 
priétés des différentielles. dQ peut s’exprimer en fonction 
de deux variables indépendantes qui peuvent étre, soit p 
et V, soit V et T, soit p et T, et on peut écrire: 

a) dQ =Adv + mwdp, 

b) dQ =Idv+ c, dT, 

c) dQ=hdp + cy dT, 
ott A, w, 1, h, cy, Cp sont les dérivées partielles de Q par 
rapport aux variables correspondantes supposées indépen- 
dantes. Mais dp, dv, dT ne sont pas indépendants, et sont 
reliés ici par 

d) pdv+ Vdp=RdT. 
On obtient alors 2 et 4 en remplagant successivement dv 
et dp par leur valeur tirée de d) et en identifiant a) et b), 
et a) et c). 

En toute rigueur, on he peut pas choisir ici T comme 
variable, puisque le gaz n’est pas en équilibre thermique et 
que T varie d’un point a l’autre de la cavité. Mais les 
expressions de 2 et y restent toujours valables. 
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Si nous nous placons dans le cas ou l’énergie 
apportée est une fonction cissoidale du temps, 


dE = E, elt dt 


la pression p varie de maniére cissoidale autour 
d’une valeur moyenne 


P= Pot piel, 
dp = io p, ei dt 


et les formules précédentes (2) et (3) s’écrivent 
respectivement: 


pix=y E\fioV, (2’) 
Pi = (y—1) Es fioV. (3’) 


La comparaison des formules (2’) et (3’) 
montre une analogie tout a fait frappante entre 
les valeurs de la pression produite respectivement 
par le thermophone et par le pistonphone. Elle 
montre aussi qu’en premiére approximation la 
pression produite par le thermophone ne dépend 
en aucune maniére des caractéristiques thermi- 
ques des feuilles d’or ou du gaz (a part le rapport 
7 = Cp/ey). 

Toutefois, si on considére les formules (2’) et 
(3’), on constate qu’elles différent par la valeur 
d’un «coefficient de rendement» en quelque sorte, 
égal a y dans un cas, et a y — 1 dans l’autre. II 
semble donc que, dans le cas du thermophone, il 
y ait une perte de l’énergie utilisable. Cela est 
bien conforme aux principes de la thermodyna- 
mique, et nous préciserons ultérieurement com- 
ment s’effectue cette perte d’énergie. 


2. Etude du mécanisme de la production de la 
pression acoustique 


La théorie précédente nous a permis de déter- 
miner rapidement la pression produite a l’inté- 
rieur de la cavité du thermophone, en appliquant 
seulement des principes généraux. Mais ce calcul 
ne donne aucune indication sur l’aspect physique 
des phénoménes. Nous allons alors examiner plus 
en détail ce qui se passe a l’intérieur de la cavité. 

Remarquons que, jusqu’a présent, nous n’avons 
émis aucune hypothése sur la maniére dont s’ef- 
fectuent les échanges de chaleur entre les feuilles 
d’or et le gaz, ou entre les différentes particules 
de gaz. Nous allons supposer maintenant que ces 
échanges se font uniquement par conduction, ce 
qui est le cas pour le domaine de température 
rencontré. Pour simplifier le probléme, nous 
allons supposer de plus qwils s’effectuent unique- 
ment par tranches paralléles au plan des feuilles 
d’or. La température du gaz en un point sera 
done fonction uniquement du temps et de la dis- 
tance x de ce point aux feuilles. 
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Nous allons considérer une particule limitée 
par deux surfaces égales et paralléles aux feuilles, 
de grandeur s, d’abscisses respectives x, et 
x + Ax. Nous prendrons ici comme variables, non 
plus V et p, mais T et p 3. Lorsqu’on passe d’un 
état défini par T et p a l'état voisin défini par 
T+ dT, p+ dp, la quantité de chaleur fournie 
au systéme est donnée par 


dQ = cpdT + hdp, 


ty (cp 9) Pe ~ 
avec h = (op/etye = R a= i 


Mais, d’aprés la théorie de la conduction, la 
quantité de chaleur dQ est égale a 


yd oT 
Os a+ Ks(7) dt 
e ale Ox |e + vx 
oT 
ce qui donne 
o2T V 


Si nous nous placgons dans le cas du régime 
cissoidal permanent, la température T(x) varie 
cissoidalement autour d’une valeur moyenne 


T= To Tete 


ou T; ne dépend que de x, et léquation (4) 
devient 


d2T, __ pia {ie AY 1@ pi 


- dx2 V Aj 


Pour simplifier l’écriture, nous poserons ~? = 
cpw/2V K. L’équation précédente s’écrit alors: 


1 EN 2 To. ee 

2ia® dx? y Po 
La solution de cette équation nous donne la 
répartition de la température en fonction de la 
distance x aux feuilles. Cette solution est valable 
pour un point situé en regard des feuilles et prés 
des feuilles. Nous supposerons que les parois de 
la cavité sont trop éloignées pour modifier sensi- 
blement les phénoménes dans cette zone. La solu- 
tion de éq.(4’) qui répond alors aux conditions aux 
limites, 4 savoir que T, doit étre égal sur les 
feuilles d’or a la température 6; de ces feuilles, 

s écrit: 

Pes (0. — yah Ty 1) oR agiay pil Ti 
rr Po Po 


* Ce qui est valable ici, puisque la particule est en équi- 
libre thermique. 
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To ew est petit de- 
Po 
vant 4), et nous le négligerons. 


D’ou 
= —1 
meepe Vier 2 TP 65) 
Y Po 
L’interprétation de l’expression (5) va nous 
étre utile pour comprendre le comportement du 
gaz. Cette expression se compose de deux 
termes. 
Considérons d’abord le premier terme 


T= 6) e—Viiex 


Dans la pratique, 


ou, si ’on revient aux variables réelles 
— ar 
Ty’ = Oi e cos (wt — xx). 


L’échauffement périodique des feuilles provo- 
que une onde thermique (comme l’ont signalé 
d’ailleurs ARNOLD et CRANDALL [2]). Cette onde 
se propage avec la vitesse w/x, variable avec la 
fréquence, puisque « est proportionnel seulement 


a Vw. De plus elle est fortement amortie, puisque 
son amplitude devient négligeable dés que la 
distance aux feuilles d’or dépasse une longueur 
d’onde 2x/x. Comme, dans le cas envisagé, x est 
voisin de 1,4 )/f, on voit que, dés que f dépasse 
une centaine de périodes par seconde, cette 
longueur d’onde est bien inférieure a la hauteur 
de la cavité, et on peut considérer que l’onde 
thermique émanant des feuilles est localisée dans 
une petite zone V; ayant la forme d’un parallélé- 
pipéde aplati situé au voisinage des feuilles (Fig.3). 


Feuilles 


ZZ. 


La 


Fig. 3. Comportement schématique de la cavité. 
. 


Le deuxiéme terme T," = “—— To Pls une 
if Po 

interprétation physique immédiate. I] entraine 
dv/V = dp/yp. Ce terme correspond donc a une 
variation adiabatique de l’état du gaz. Ce terme 
est une constante et ne dépend donc pas de la 
position du point considéré dans la cavité. Il est 
done valable, méme si on s’éloigne des feuilles et si 
on quitte le domaine ot l’équation (4’) est valable. 
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En définitive, le phénoméne de production du 
son a l’intérieur de la cavité s’effectue schémati- 
quement de la maniére suivante: Les variations 
de température des feuilles d’or donnent nais- 
sance a une onde thermique fortement amortie et 
localisée (pour des fréquences pas trop basses) a 
Vintérieur d’une mince couche gazeuse V1 avoisi- 
nant les feuilles d’or. Sous Vinfluence de cette 
onde, les particules de gaz de V; sont soumises 
a des variations de température, donc a des varia- 
tions de volume. II s’ensuit une variation pério- 
dique de l’épaisseur de la couche et tout se passe 
comme si la surface de cette couche se déplagait 
a la maniére d’un piston aérodynamique. Sous 
Pinfluence de ce piston, la masse de gaz occupant 
le reste de la cavité, V2, est soumise a des varia- 
tions de volume, tout comme dans le cas d’un 
pistonphone ordinaire, et elle se comporte de 
maniére uniforme et adiabatique: La variation 
de volume de V2 due a un déplacement du piston 
aérodynamique entraine une variation de pres- 
sion. On distingue donc nettement dans le com- 
portement du gaz deux régions essentiellement 
différentes. L’une située au voisinage des feuilles 
d’or, ot les phénoménes thermiques sont prépon- 
dérants, et l’autre constituée par le reste de la 
cavité, ou les échanges thermiques sont nuls et ou 
se produit un phénoméne de compression adia- 
batique et uniforme: Le thermophone peut étre 
considéré comme un pistonphone aérodyna- 
mique. 

Considérons alors de nouveau deux états voi- 
sins du volume V. Ces deux états correspondent, 
pour le systéme Vj, a des pressions p et p+ dp, 
et des volumes V; et V+ dv, et, de méme, pour 
le systéme V2, a des pressions p et p+ dp, et des 
volumes V2 et V2—dv, le dp et le dv étant res- 
pectivement les mémes pour les deux systémes. 

Le systéme V2 se comporte comme un piston- 
phone ordinaire. La variation de pression dp 
correspond a un travail dr du piston aérodyna- 
mique donné, d’aprés (2), par 


V2 A 
dr = — dp s — dp. 
ie - Ie P 


Mais cette méme variation de pression dp 
correspond, pour le systéme Vy, 4 l’échange de 
Vénergie dE avec le milieu extérieur, donnée, 
d’aprés éq. (3), par 

V 
dE = ———- dp. 
pa 

La comparaison des deux formules précédentes 
conduit a. 


dr = — dk. (6) 
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Il y a donc seulement une fraction de l’énergie 
fournie par l’extérieur au systéme qui est trans- 
formée en travail. Ce n’est que cette fraction qui 
est utilisée pour produire la variation de pression 
dans la cavité. Le reste est dissipé dans le volume 
V, sous forme de chaleur. On retrouve le fait bien 
connu de la dégradation de l’énergie. 

On peut toutefois étre surpris de constater que 
la fraction de l’énergie utilisable, égale a (y—1)/y, 
soit 2/7 environ, soit constante et ne dépende ni 
de la cavité, ni des propriétés thermiques du gaz. 
De fait, on peut expliquer cela de la maniére 
suivante: La variation de pression dp correspond 
pour le systéme V2, donc pour le systéme Vj, a 
une variation de volume dv donnée par: 

dp V2_dpV_y—ldE_ RdkE_ dE 

jae ann ce me ee ne ea a8 
La comparaison de cette derniére expression 
avec éq. (1) montre que la variation de volume de 
V, est pratiquement la méme que si la transfor- 
mation du systéme Vj  s’effectuait A pression 
constante, c’est-a-dire que si V2 était infini. 

Nous avons pu décomposer dans l’espace la 
cavité en deux systémes distincts fonctionnant 
suivant des modes tout a fait différents. Cette 
décomposition n’est valable en toute rigueur que 
pour des fréquences relativement élevées. Pour 
les basses fréquences, la limite entre V; et V2 
s’estompe, et l’hypothése selon laquelle V2 est 
grand devant Vj, n’est plus valable. On ne peut 
plus distinguer une zone ou l’énergie calorifique 
se transforme en travail, ce travail se transmet- 
tant au reste de la cavité par l’intermédiaire d’un 
piston aérodynamique pour donner naissance a une 
pression acoustique. Cependant, et ceci est tout a 
fait remarquable, quelle que soit la fréquence, 
il y a toujours la méme fraction, égale a (y — 1)/y, 
de l’énergie apportée au systéme, qui est utilisée 
pour donner naissance 4 une pression acoustique. 


dv 


3. Théorie approchée 


La théorie précédente est évidemment un peu 
sommaire. Nous avons supposé en effet que toute 
V’énergie apportée par effet Joule se dissipait dans 
le gaz et que, de plus, il n’y avait pas de perte de 
chaleur par les parois de la cavité. Ce qui compli- 
que beaucoup la théorie du thermophone, c’est 
que les feuilles d’or absorbent une partie impor- 
tante de l’énergie calorifique, et que, comme la 
cavité est petite, les effets de surface ne sont pas 
négligeables. Pour le pistonphone, au contraire, 
la formule (2’) s’applique pratiquement sans 
correction. Nous allons examiner alors successive- 
ment l’influence des feuilles d’or et celle des parois 
de la cavité. 
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a) Effet des feuilles d’or 


Si nous considérons le systéme formé par le 
volume gazeux V, la théorie du premier para- 
graphe reste entiérement applicable, et les for- 
mules (3) et (3’) sont toujours valables, a condi- 
tion de considérer que les quantités dE et Ej; 
correspondent a l’énergie absorbée par le gaz, 
mais non plus a l’énergie fournie par effet Joule. 
Il faut donc chercher a calculer l'une en fonction 
de l’autre. 

Considérons toujours deux états voisins. Lors- 
que le systéme passe.de l'état 1 a état 2, le 
milieu extérieur a apporté par effet Joule une 
quantité d’énergie dW, mais les feuilles se sont 
échauffées, et leur température s’est accrue de 
dé, ce qui a nécessité une quantité de chaleur 
ay d6, ou a est la surface des feuilles, et y leur 
capacité calorifique par unité de surface. On a 
donc 


dE = dW— Jay do. (7) 


Mais, si on suppose que l’échange de chaleur 
entre le gaz et les feuilles se fait par conduction, 


x 


aE = 25 Ka(Z) dt. (8) 
x=0 


Si nous nous plagons dans le cas du régime 
cissoidal permanent, (7) est remplacé par 


E, = W,—Jaionrzhy 


et (8) est de méme remplacé par 


(7’) 


(8’) 
Mais T}, solution de (4’) est donnée par (5). D’ot 


(dT, se 
(a Wher = V2i« 0, 
et Fy =2J Kal/2ix0)=2JKa(1+i)aby. (8”) 


En éliminant 6; entre (7’) et (8”), on trouve 
finalement 

Wy 
(l+i)zyo ° 

4K x 


Ey (9) 


Et 


b) Effet des parois 


D’aprés les considérations de la deuxiéme par- 
tie, effet des parois n’est sensible que pour le 
volume V2, et on peut traiter le probléme comme 
dans le cas d’un pistonphone ordinaire. Considé- 
rons alors une cavité de volume V soumise a des 
compressions adiabatiques périodiques, et sup- 
posons que les parois de la cavité soient conduc- 
trices de la chaleur. A V’intérieur de la cavité, les 
variations de volume du gaz entrainent des varia- 
tions de température, mais, sur les parois, la 


———eer hr 
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température ne varie pas. Il y a donc une diffé- 
rence périodique entre la température de la 
masse du gaz et celle des parois. Il va donc se 
passer un phénoméne tout a fait analogue a celui 
qui se passe au voisinage des feuilles, c’est-a- 


dire la naissance d’une onde thermique qui va ° 


provoquer des variations de volume d’une couche 
gazeuse située au voisinage des parois de la cavité, 
et ces variations de volume modifient la pression 
acoustique. 

Si 6,’ désigne la variation de température du 
gaz au milieu de la cavité, la variation de tempé- 
rature a une distance x; des parois sera donnée 
sensiblement par 


Oy. (1 —e-— V2ian), 


Or la formule (8”) nous donne une relation 
entre l’énergie FE, apportée au systéme, et la 
température 0; du gaz contre les feuilles d’or. 
On en déduit par analogie que l’énergie absorbée 
par les parois de la cavité, de surface S, est égale a 


Ey’ =JKS/2i« 6’. 


Mais 6;’ est la variation de température due a 
une transformation adiabatique donnant nais- 
sance a une pression acoustique p1 


—] V 
haf aig je ssl 
Y Po J cp 
L’énergie réelle fournie au systéme est donc 
finalement égale a 
Wy 


(1+i)xyo 
ba 4K « 


Ek, = Ey’. 


Quant a la pression produite, elle est toujours 


donnée par (3’) 


THERMOPHONE 


Belice d jy 
area aaa 
d’ot 
, (y—1) My, y—1 KS(1+i)apiV 
1=>- = . 
’ (l+i)yo| ioV Cy 
1@ y[1+ eee 
ce qui donne finalement (10) 
vipa (y—1) Mi 
pig (l+i)yo esi 1) 
io VL 4K x ad 2Va 


L’expression (10) donne la valeur réelle de la 
pression produite. On retrouve la formule dé- 
montrée d’une maniére plus rigoureuse dans 
Varticle paru dans les Annales des Télécommuni- 
cations. 
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En résumé, le fonctionnement du thermophone 
se raméne au schéma suivant (Fig. 4): Une source 


festa faite, 


Générateur 


Fig. 4. Fonctionnement schématique du thermophone. 


d’énergie extérieure (électrique) fournit une 
quantité d’énergie dW. Une partie de cette 
énergie (J dQ’) est perdue et absorbée pour 
échauffer les feuilles d’or. Seule la quantité (J dQ) 
est transmise a la cavité. Une partie de cette 
énergie calorifique est perdue par effet thermique 
dans le volume V; du gaz, le reste est communi- 
qué sous forme de travail au piston P, dans la 
proportion 
dt vil 
of SK Cred heey 

L’énergie sous forme de travail dt est alors 
transmise au volume V, et sert a augmenter 
l’énergie interne du gaz (augmentation de tempé- 
rature et augmentation de pression). Seule une 
petite fraction (J dQ") est dissipée sous forme de 
chaleur a travers les parois de la cavité. Au total 
on a une transformation d’énergie électrique en 


énergie calorifique, puis d’énergie calorifique en 
travail, puis de travail en énergie interne. Toutes 
ces transformations s’effectuent avec perte (dans 
les feuilles, dans le volume V,, dans les parois) 
et finalement le rendement est médiocre, surtout 
dans les hautes fréquences, c’est-a-dire lorsque 
la transformation s’effectue rapidement, car, a 
ce moment 1a, une partie importante de l’énergie 
est perdue dans les feuilles d’or, les calories 
n’ayant pas le temps de se dissiper par conduc- 
tion dans le volume Vj. 


Conclusion 


Nous avons essayé, au cours de cet exposé, de 
dégager les aspects essentiels du fonctionnement 
du thermophone. Grace a des considérations 
presque purement physiques, nous sommes arri- 
vés 4 calculer la pression acoustique développée 
dans la cavité. La formule a laquelle nous avons 
abouti est identique a-la formule trouvée d’autre 
part a l’aide de considérations théoriques plus 
rigoureuses. Le schéma que nous avons tracé est 
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donc suffisamment valable, bien que nous ayons 
négligé certains phénoménes qui s’avérent finale- 
ment de peu d’importance. 

Si nous nous plagons maintenant d’un point 
de vue pratique, a savoir, de l'utilisation du 
thermophone pour ]’étalonnage des microphones, 
quelle conclusion allons-nous pouvoir tirer ? On 
sait que la méthode du thermophone a eu de 
nombreux détracteurs. Il lui a été reproché 
notamment la complexité et linsuffisance de sa 
théorie et l’incertitude sur les propriétés thermi- 
ques des feuilles d’or et du gaz. Il semble que, 
sans vouloir réhabiliter tout a fait le thermophone, 
qui présente sans conteste des inconvénients, 
l’exposé que nous venons de faire apporte quel- 
ques éclaircissements sur ce point. D’une part, 
et sans prétendre avoir épuisé le probléme, nous 
pensons pouvoir affirmer que la théorie exposée 
dans les Annales des Télécommunications a été 
établie avec assez de rigueur et de précision pour 
étre admise comme physiquement valable, au 
moins pour des fréquences pas trop basses. 
D’autre part, comme les feuilles d’or apparaissent 
dans cette théorie comme un élément secondaire 
et correctif, les incertitudes sur leurs propriétés, 
thermiques ont également une influence secon- 
daire dans le calcul de la pression. De fait il 
s’avére dans la pratique que le thermophone reste 
un instrument d’étalonnage utilisable et qui a 
déja rendu de nombreux services. Toutefois cette 
influence se fait sentir malgré tout, et de maniére 
sensible, dés que la fréquence est un peu élevée, 
et elle limite par suite les possibilités de l’appa- 
reil. Nous ne pensons pas que dans les conditions 
actuelles, la précision obtenue avec le thermo- 
phone puisse étre meilleure que 0,5 dB, méme en 
utilisant plusieurs thermophones pour faire la 
mesure, alors que la méthode de réciprocité 
permet d’atteindre en principe une précision de 
0,2 dB; et il n’est pas douteux que cette derniére 
méthode, quoique longue, ne soit utilisée de plus 
en plus largement. 

Le travail que nous avons entrepris présente 
done un intérét plutét théorique que pratique. 
I] avait surtout pour but de tenter d’expliquer 
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les écarts trouvés par différents auteurs [7], [8], 
entre les résultats obtenus par la méthode du 
thermophone et d’autres méthodes classiques 
(pistonphone, réciprocité, disque de Rayleigh). 
Bien que nous n’ayons pu faire encore une expéri- 
mentation suffisamment précise et compléte pour 
pouvoir tirer des conclusions définitives, nous 
avons déja obtenu des résultats assez encoura- 
geants pour pouvoir affirmer que la nouvelle 
théorie du thermophone se vérifie bien expéri- 
mentalement et qu’on trouve un accord accep- 
table (compte-tenu des erreurs de mesure), entre 
les différentes méthodes. 

Ce travail nous a permis en outre de mettre en 
évidence le fonctionnement du thermophone du 
point de vue physique, et il nous a paru intéres- 
sant de l’exposer ici en détail. Nous avons vu que 
ce fonctionnement d’apparence complexe pouvait 
s’expliquer simplement par une succession de 
transformations énergétiques élémentaires, et il 
est tout a fait remarquable que ces simples consi- 
dérations énergétiques nous aient permis de cal- 
culer de maniére suffisamment approchée la 


pression acoustique produite dans la cavité. 
(Regu le 31 Octobre 1955.) 
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UBER DEN STROMUNGSWIDERSTAND POROSER BAUSTEINE* 


von F. BARTHEL, Berlin 


Zusammenfassung 


Es werden der statische Strémungswiderstand und die Porositat an den gebrauchlichsten 
porésen Bausteinen gemessen. Man findet ein zum Teil sehr starkes Anwachsen des Strémungs- 
widerstandes mit der Strémungsgeschwindigkeit, wahrend die untersuchten feinporigen Schall- 
schluckstoffe diese Eigenschaft nicht aufweisen. Beim Vergleich des so gemessenen Strémungs- 
widerstandes mit dem auf Grund von Wandwiderstandsmessungen im Rohr sich ergebenden 
akustischen Strémungswiderstand erhalt man bei den feinporigen Schallschluckstoffen gute 
Ubereinstimmung, dagegen liegt bei Proben aus pordsen Bausteinen der akustische Strémungs- 
widerstand stets hoher als der statische und wichst etwa oberhalb eines Schallpegels von 100 dB 
mit der Schnelle. Bei den untersuchten Bausteinproben werden auferdem Werte des Struk- 
turfaktors bis zu 10 gemessen. 


Summary 


The static flow resistance and porosity of common porous building materials are measured. 
The flow resistance in some cases rises very steeply with air velocity, a result not found for 
the measured sound absorbing materials with fine pores. By comparing the measured resistance 
with acoustic resistance, deduced from measurements of wall impedances in a tube, one gets good 
agreement for absorbent materials with fine pores, but in samples of brick acoustic flow resistance 
is larger than static flow resistance and rises beyond a sound pressure level of 100 dB with air 
velocity. Values of “‘structure factor” up to 10 are found for the brick samples. 


Sommaire 


On a mesuré la résistance statique 4 l’écoulement et la porosité des matériaux poreux de cons- 
truction les plus courants. On a trouvé que cette résistance augmente souvent fortement avec la 
vitesse d’écoulement, résultat que ne donnent pas les matériaux absorbants mesuré: a petits pores. 
Si on compare la résistance a l’écoulement ainsi mesurée 4 la résistance acoustique a |’écoulement 
qui est obtenu par des mesures de résistance de paroi faites dans un tube, on constate un accord 
satisfaisant dans le cas des matériaux absorbants 4 petits pores; dans le cas des matériaux de 
construction poreux, par contre, la résistance acoustique 4 "écoulement est toujours supérieure 


a la résistance statique, et augmente a plus d’environ 100 dB quand l’amplitude de la célérité 
croit. On a trouvé, en outre, pour les matériaux de construction qui ont été étudiés, des valeurs 


du facteur de structure allant jusqu’a 10. 


1. Einleitung 


Uber die Messung des Strémungswiderstandes 
von porésen Schallschluckstoffen ist eine ausge- 
dehnte Literatur vorhanden [1]---[7]. In jiingster 
Zeit wurde vom Deutschen Fachnormenaus- 
schuB Materialpriifung nach Vorschlagen von 
G. VenzKeE ein Normungsentwurf aufgestellt [8]. 
Der sogenannte auBere Strémungswiderstand W 
einer Probe der Dicke d ist definiert als das Ver- 
haltnis von der Druckdifferenz Ap zur Stré- 
mungsgeschwindigkeit v vor und hinter der 
Probe, wenn durch sie ein konstanter, nur in 
einer Richtung erfolgender Luftstrom hindurch- 
flieBt. Es ist also W= 5d = Ap/v, wobei £ der 
langenspezifische Strémungswiderstand (Stré- 
mungskennwiderstand) ist. 

Befindet sich eine Materialprobe mit der Poro- 
sitat o und dem Strukturfaktor 7 in der Entfer- 


* Auszug aus einer Diplomarbeit des I. Physikalischen 
Instituts der Universitat Miinchen. 


nung / vor dem starren Abschlu8 eines Kundtschen 
Rohres, dann gilt fiir den Wandwiderstand W vor 
der Probe, falls die Probendicke d sehr klein 


gegen die Wellenlange des benutzten Tones ist, 
WE=f£d+ Eva 
o 


dabei ist 9 die Luftdichte, Z der Kennwiderstand 
der Luft und k= w/c, wobei w= 27f die Kreis- 
frequenz des benutzten Tones und c die Schall- 
geschwindigkeit in Luft ist. Man kann also aus 
einer Wandwiderstandsmessung den Strémungs- 
widerstand bestimmen, und zwar ist dieser im 
Resonanzfall, d. h. bei verschwindendem Imagi- 
narteil, gleich dem Wandwiderstand selbst. Es 
ist bekannt, daB im allgemeinen bei den Schall- 
schluckstoffen der akustische und statische Stré- 
mungswiderstand gut iibereinstimmen. Bei grob- 
porigen Substanzen, wie etwa bei porésen Bau- 
steinen, sind aber Unterschiede zwischen beiden 
Werten zu erwarten. 
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2. Statische Messungen an porosen Bausteinen 


a) MefBapparatur 

Es wurde mit einem elektrisch betriebenen 
Geblase ein kontinuierlicher Luftstrom von 
0,5 bis 20 1/min bei einem Gegendruck von héch- 
stens 20 cm Wassersaule erzeugt. Der Luftstrom 
durchlief eine Beruhigungsstrecke und wurde mit 
einem grofBen und einem kleinen Steuerventil 
grob- bzw. feingeregelt und dann durch den 
DurchfluBmesser, einen nassen Kontrollgaszah- 
ler, geleitet (Abb. 1). Nach dem DurchfluBmesser 
strémte die Luft in ein Metallrohr von 40 cm 
Lange und 86 mm Durchmesser ein und wurde 


Filz 


Gaszdhler 


Geblase 


Abb. 1. Schema der Apparatur zur Messung des statischen 
Strémungswiderstandes. 


dabei durch eine 2 cm dicke Filzschicht in eine 
laminare Strémung verwandelt. An dieses Rohr 
wurde ein Metallmantel von quadratischem 
Querschnitt (15 x 15 cm?) und 20cm Mantel- 
lange angeschlossen, in den man die Baustein- 
proben eingipste. Der Ubergang vom kreisfér- 
migen zum quadratischen Querschnitt erfolgte 
dabei mittels einer Metallplatte vom Querschnitt 
20 x 20 cm2, welche in der Mitte eine kreisrunde 
Offnung mit dem Durchmesser des Rohres 
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(86 mm) hatte, an das sie angeflanscht wurde. An 
diese Metallplatte konnte der Halterungsmantel 
mit der Bausteinprobe dann direkt angeflanscht 
werden. Die Bausteinproben hatten den Quer- 
schnitt 11 x 11 cm? und eine Starke von 2 bis 5 em. 
Sie waren aus beliebigen Steinen einer Fabrika- 
tionsserie mit einer Spezialkreissige, der soge- 
nannten Karborundsiage, herausgeschnitten wor- 
den und wurden in den Metallmantel so ein- 
gegipst, dafi méglichst wenig feuchter Gips in 
das Innere des Steines eindringen konnte. 

Es wurde dann mit einem Mikromanometer 
das Druckgefalle gemessen, das den konstanten 
Luftstrom, dessen Geschwindigkeit aus der An- 
zeige des DurchfluBmessers hervorging, durch die 
untersuchte Probe trieb. Das Mikromanometer 
und der DurchfluBmesser wurden geeicht. Beim 
Mikromanometer, das nach dem Prinzip des 
Schragrohr-Fliissigkeitsmanometers arbeitete, er- 
folgte die Eichung mittels einer Druckwaage. Es 
ergab sich eine MeBunsicherheit des Mikromano- 
meters von 0,01 mm Wassersdule. Bei der 
Eichung des Gasmessers lieB man aus einem 
KubiziergefaB mittels elektrischer Kontakte 
genau 1001 Luft durch den Gasmesser strémen 
und beobachtete dabei, wieviel Liter von dem 
Gasmesser angezeigt wurden. Zur Bestimmung 
der Porositét der untersuchten Proben konnte 
man mit einer Feinwaage den Auftrieb der Pro- 
ben in Alkohol messen und daraus die Porositat 
errechnen. 


b) MefBergebnisse 

Die Messungen wurden an Bims-, Tuff-, 
Schlacken-, Ziegelsplitt- und Schwerbetonsteinen 
vorgenommen, welche die gebrauchlichsten poré- 
sen Bausteinarten sind. Bei Tuff- und Schlacken- 
steinen wurden Proben aus je zwei Serien gemes- 


Tabelle I 


Materialart 


Bims econ te iets 0,37 
Bims Wires eee 0,39 
Tuff, Serie Tae. 0,365 
Taft Semmes: T tise «ct 0,39 
Teil. ened Esc. 5 0,335 
‘Toft. Serietl Lob... 0,315 
Schlackenstein I ... 0,29 
Schlackenstein II .. 0,25 
Schlackenstein II .. 0,235 
Ziegelsplitt ....... 0,3 
Ziegelsplitt........ os 
Schwerbeton ...... 0,27 
Schwerbeton ...... 0,16 


statischer spezifischer Str6mungs- 
widerstand = in g em~$ st 


bei 1,5 cm/s 


bei 0,3 cm/s bei 3,0 cm/s. 


P 
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sen, um den Unterschied der verschiedenen Her- 
stellungsserien bei gleicher Steinart festzustellen. 
In der Tabelle I findet man die Werte der Poro- 
sitat, des Raumgewichtes sowie des spezifischen 
statischen Strémungswiderstandes bei den Stré- 


mungsgeschwindigkeiten 0,3, 1,5 und 3 cm/s. | 


In Tabelle I sind zu jeder Serie die MeBwerte von 
jeweils zwei Bausteinproben angegeben, die aus 
den Steinen der Serie herausgesigt wurden. Bei 
Schlackenstein Serie I war nur eine Probe ver- 
fiigbar. Mit Ausnahme des Ziegelsplittsteines 
waren alle Steine Wandsteine, wahrend der Zie- 
gelsplittstein ein Deckenfillkérper war und des- 
halb eine bedeutend dichtere Bindeschicht hatte, 
da er fester und fetter abgebunden sein muBte. 
Das Zuschlagmaterial hatte bei den Steinen der 
Schatzung nach eine mittlere Korngrésse von 
4 bis 8mm Korndurchmesser. Im Ziegelsplittdek- 


3 cmis 
Strémungsgeschwindigkeit ——= 


Abb. 2. Spezifischer Str6mungswiderstand von zwei Bims- 
steinen aus gleicher Herstellungsserie. Aus jedem 
Stein wurden zwei Proben gemessen; 
x x 1. Stein, @ --- @ 2. Stein. 


2000 


gcm-3s- 


1500 


500 


0 1 2 cmis 3 
Stromungsgeschwindigkeit ——=— 


Abb. 3. Spezifischer Str6émungswiderstand von zwei Zie- 
gelsplittdeckensteinen aus gleicher Herstellungs- 
serie. Aus jedem Stein wurden zwei Proben ge- 
messen, x x 1. Stein, @ --- @ 2. Stein. 
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kenstein war auBerdem noch ein etwas gréberer 
Zuschlag zur Erhéhung der mechanischen Festig- 
keit vorhanden. Zur Erginzung zu den in der 
Tabelle I festgehaltenen Messungen wurden noch 
bei Bimsstein und Ziegelsplittstein aus zwei 
Steinen der untersuchten Serie jeweils zwei Pro- 
ben gemessen, um zu sehen, wie grof} die Streu- 
ung fiir Proben ist, die aus dem gleichen Stein 
sind, und um die Streuung der einzelnen Steine 
einer Serie untereinander zu untersuchen (Abb. 2 


und Abb. 3). 


c) Diskussion der Messungen 


Der spezifische Strémungswiderstand der poré- 
sen Bausteine wird hauptsachlich durch die 
Bindeschicht bestimmt, welche die Kérner des 
Zuschlagmaterials untereinander verkittet. Dem- 
gegeniiber ist die Art und KorngréBe des Zu- 
schlagmaterials von geringerer Bedeutung. Ist 
die Bindeschicht porig, so erhalt man den ver- 
haltnismaBig geringen Strémungswiderstand der 
Tuff-, Bims- und Schlackensteine. Macht man 
aber bei der Fabrikation der Steine die Binde- 
schicht durch nasse und fette Abbindung oder 
durch Riitteln, Pressen usw. dichter, so erhéht 
sich der Strémungswiderstand erheblich, wie das 
bei Schwerbeton und besonders bei Ziegelsplitt 
der Fall war. 

Fur die gemessenen Tuffstein-, Bimsstein- und 
Schlackensteinsorten laBt sich zusammenfassend 
sagen, daf} der statische Strémungswiderstand 
bei Dicken von wenigen Zentimetern, wenn man 
den fiir die Geschwindigkeit Null extrapolierten 
Wert nimmt, nur wenig vom Kennwiderstand 
der Luft verschieden war. Daraus folgt, daB 
diese Materialien, was ihren statischen Stré- 
mungswiderstand angeht, recht gut zur Schall- 
schluckung geeignet sind. Allerdings war ihre 
Porositat mit Werten zwischen 0,4 bis 0,25 wesent- 
lich kleiner als bei den gewéhnlichen Schall- 
schluckstoffen, deren Porositat etwa bei 0,8 liegt. 
Bei diesen Bausteinen blieb auch die Material- 
streuung der einzelnen Proben innerhalb einer 
untersuchten Fabrikationsserie in tragbaren Gren- 
zen, so das durch Untersuchung mehrerer Pro- 
ben ein Mittelwert fiir den spezifischen Stré- 
mungswiderstand der Serie erhalten werden 
kann, fiir den die Streuungen der Kinzelmessung 
héchstens + 50% betragen. Dagegen war bei den 
Ziegelsplittdeckensteinen die Materialstreuung 
so stark, da} man offenbar zu gréBeren Proben- 
querschnitten greifen miBte, wenn man den 
spezifischen Strémungswiderstand als Material- 
konstante und nicht, wie hier, zum physikali- 
schen Uberblick ermitteln méchte. 
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Wie sich aus Tabelle I ergibt, war die Material- 
streuung innerhalb einer Sorte bei den Messungen 
der Porositat und auch des Raumgewichtes ver- 
haltnismaBig gering, jedenfalls geringer als bei 
den Strémungswiderstandsmessungen. Erstaun- 
lich ist, daB mit Ausnahme der einen Schwer- 
betonprobe die Porositat aller gemessenen Stein- 
proben in den engen Bereich von 0,23 bis 0,39 fiel. 
Die Tatsache, da die Porositat bei den Tuff-, 
Bims- und Schlackensteinen etwa die gleiche 
war wie bei den Ziegelsplittdeckensteinen, obwohl 
letztere eine wesentlich dichtere Bindeschicht 
hatten, erklart sich daraus, daf} bei ersteren 
zwar viele, dafiir aber kleine Hohlraume vor- 
handen waren, wahrend die Ziegelsplittsteine 
weniger aber dafiir gréBere Hohlraume hatten. 

Es wurde bei allen gemessenen Bausteinen 
infolge der starken Querschnittsveranderungen, 
welche die Luft in den Kandlen und Hohlraumen 
der Bausteine zu durchstrémen hatte, ein mehr 
oder minder starkes Anwachsen des Strémungs- 
widerstandes mit der Geschwindigkeit fest- 
gestellt. Vermutlich bildeten sich an den Quer- 
schnittserweiterungen Wirbel, welche die Nicht- 
linearitat der Abhangigkeit der Druckdifferenz 
von der Strémungsgeschwindigkeit zur Folge 
haben. Bemerkenswert dabei ist, daB bei zahl- 
reichen Proben eine lineare Abhangigkeit des 
Strémungswiderstandes von der Geschwindigkeit 
in dem untersuchten Bereich von 0,3 bis 3 cm/s 
festgestellt wurde. 

Besonders stark trat die Nichtlinearitat bei der 
Tuffsteinserie I, bei der Schlackensteinserie I, bei 
dem Schwerbeton und bei dem Ziegelsplitt- 
deckenstein in Erscheinung. Die Steine mit dem 
gréBeren spezifischen Strémungswiderstand wie- 
sen im allgemeinen einen starkeren Anstieg des 
Strémungswiderstandes mit der Geschwindigkeit 
auf, vermutlich weil sie weniger durchgehende 
Kaniale hatten, in denen daher die Strémungs- 
geschwindigkeit betrachtlich ttberhéht wurde. 


3. Statische und akustische Vergleichsmessungen 
a) Mefapparatur 


Bei den Vergleichsmessungen befanden sich 
die untersuchten Proben in einem Halterungs- 
rohr mit 85 mm Innendurchmesser und 20 cm 
Rohrlange. Das Halterungsrohr konnte sowohl 
an das Strémungsrohr der im vorigen Abschnitt 
beschriebenen Apparatur zur statischen Bestim- 
mung des Strémungswiderstandes wie an ein 
Kundtsches Rohr angeschlossen werden, wobei 
das Halterungsrohr auch den gleichen Innen- 
durchmesser wie diese beiden Rohre hatte. Man 
konnte so statisch und akustisch den Strémungs- 
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widerstand der Proben messen, ohne ihre Rand- 
dichtung zu verandern. Bei den untersuchten 
Glaswolleschichten wurden die Proben aus den 
Glaswollematten ausgeschnitten und sorgfaltig 
in das Halterungsrohr eingepabt. Bei den porésen 
Bausteinen wurden die Proben in den Fabriken 
aus der Mischung genommen, die zur Serienher- 
stellung diente und in das Halterungsrohr ein- 
gefullt. Nach der vorgeschriebenen Trockenzeit 
konnten dann die statischen und akustischen 
Messungen vorgenommen werden. Die akustische 
Bestimmung des Strémungswiderstandes erfolgte 
nach der in der Finleitung angegebenen Formel 
durch die Messung des Wandwiderstandes vor der 
Probe, wenn diese mit Luftabstand vor dem star- 
ren Abschluf angebracht war. Der gewiinschte 
Abstand der Probe von dem starren Rohr- 
abschluB wurde dadurch erreicht, dafs man eine 
1 cm starke Metallplatte mit eingelegter Gummi- 
dichtung in dem Halterungsrohr, das an das 
Kundtsche Rohr angeschlossen wurde, verschie- 
ben konnte. Zur Messung des Wandwiderstandes 
diente in bekannter Weise die Methode der ste- 
henden Wellen im Kundtschen Rohr. Die Zuver- 
lassigkeit des benutzten Rohres wurde dabei 
durch Ausmessen der stehenden Wellen bei star- 
rem RohrabschluB geprift, wobei insbesondere 
die durch die Rohrdimpfung sich ergebende 
Enveloppe der Schalldruckminima untersucht 


wurde [9]. 
b) Messungen an feinporigen Schallschluckstoffen 


Je eine 2 cm starke Schicht aus Sillan-Wolle, 
Gerrix-Glaswolle und  Schallex-Glasfaserstoff 
wurde in 10 em Abstand vor dem starren Rohr- 


statisch bestimmter Strémungswiderstand 


Rie 


Abb. 4. Wandwiderstandsdiagramm fiir drei 2 em starke 
feinporige Schluckstoffproben bei 10 em Abstand 
vom starren Rohrende. X ist der Imaginarteil, 
R der Realteil des Wandwiderstandes, Z der Kenn- 
widerstand der Luft. 

@—e—e Schallex-Platte, o—: 
x—x—»x Gerrix-Matte. 


»>—o Sillan-Wolle, 


Ee 
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abschluB auf ihren Wandwiderstand bei den 
Frequenzen 400, 500, 640, 800 und 1000 Hz unter- 
sucht. In Abb. 4 sind die gemessenen Wand- 
widerstande festgehalten. Man erkennt aus 


Abb. 4, da die einzelnen Wandwiderstands- 


kurven an den Stellen durch die reelle Achse’ 


gehen, an denen der Realteil gleich dem statisch 
bestimmten Strémungswiderstand ist. Damit 
herrschte bei den drei Proben innerhalb der MeB- 
streuung, die man bei den Wandwiderstands- 
bestimmungen mit 5% annehmen durfte, Uber- 
einstimmung zwischen dem statisch gemessenen 
und dem akustisch bestimmten Strémungswider- 
stand. Es erwies sich dabei sowohl der statisch 
wie akustisch gemessene Strémungswiderstand 
als unabhangig von der DurchfluBgeschwindig- 
keit, die bis zu 5 cm/s betrug, bzw. der Schall- 
intensitat, welche bis zu 130 dB gesteigert wurde. 


c) Messungen an Proben poréser Bauelemente 
a) Strukturfaktor 


Im Innern eines porésen Stoffes mit dem Struk- 
turfaktor y ist die Geschwindigkeit des Luftschal- 
les um den Faktor 1/)/7 kleiner als im freien 
Schallfeld. Da bei einem porésen Material mit 
nicht zu hohem Strémungswiderstand, das un- 
mittelbar vor starrer Wand angebracht ist, der 
Schluckgrad ein Maximum hat, wenn die Schicht- 
dicke gleich einer Viertelwellenlange wird, kann 
man aus der Frequenz f des Schluckgradmaxi- 
mums und aus der Schichtdicke d nach der For- 
mel //y = c/4.d f den Strukturfaktor bestimmen, 
c ist die Schallgeschwindigkeit im freien Luft- 
schallfeld. Der nach dieser Formel berechnete 
Strukturfaktor ist in der Tabelle II fir fiinf ver- 
schiedene Proben aus porésen Baustoffen ange- 
geben, welche alle unmittelbar vor starrem Rohr- 
abschluB auf ihren Schluckgrad untersucht wur- 
den. In Abb. 5 ist der Schluckgradverlauf in 
seiner Abhangigkeit von der Frequenz fiir drei 
von den gemessenen Proben aufgetragen. 


0 400 800 


f —e 


1600 Hz 3200 


Abb. 5. Frequenzabhingigkeit des Schluckgrades ao fiir 
drei Baustoffproben der Dicke d unmittelbar vor 
starrem RohrabschluB; 

x x Heraklithplatte, d= 2,5¢ 
e——e Schlackenstein, d= 3,3¢ 
Oo -- © Tuffstein, d= 2,2 ¢ 
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Tabelle II 


gem? s-1 


6) Messung des Strémungswiderstandes 
im geschwindigkeitsunabhangigen 
(linearen) Bereich 


Es wurden an fiinf Baustoffproben Wand- 
widerstandsmessungen und statische Strémungs- 
widerstandsmessungen vorgenommen. Die Wand- 
widerstandsmessungen erfolgten etwa im Fre- 
quenzbereich von 200 bis 800 Hz. Aus der Tat- 
sache, daf} fiir diese Baustoffproben der Reso- 
nanzfall, d.h. der rein reelle Wandwiderstand, 
schon unterhalb der Frequenzen eintrat, fir 
welche der Abstand zwischen dem starren Rohr- 
abschlu8 und den Proben gleich einer Viertel- 
wellenlange gewesen ware, kann man folgern, 
daB die Baustoffproben teilweise wie Lochplatten 
wirkten. Man kann dieses Verhalten der porésen 
Baustoffe angenahert berechnen, wenn man ne- 
ben der Schichtdicke d und dem Wandabstand I 
noch den Strukturfaktor sowie die Porositat der 
Proben kennt. Denn unter der Bedingung, dab 
der Wandabstand /| klein gegen die Wellenlange 
ist, gilt fir den Wandwiderstand W vor der 
Probe: 


a Ae. Ah SARS. a 
[} ee inzo ia t+ 5d, 

dabei ist K = 1,4-10° g cm-1s~2 die Raumsteife 
und 9 = 1,2-10-3 gcm~* die Dichte der Luft. 
In unserem Falle war z. B. fiir die 2,5 cm starke 
Heraklithprobe bei 1=6cm, y= 7,3 (nach 
obigem Verfahren gemessen) und o = 0,6 die 
berechnete Resonanzfrequenz gleich 390 Hz. Die 
gemessene Resonanzfrequenz lag allerdings bei 
280 Hz, was wohl davon herriihren diirfte, da 
die Probendicke nicht mehr als geniigend klein 
gegeniiber der durch den Strukturfaktor verkirz- 
ten Wellenlange anzusehen war. AuBerdem 
konnte die Porositat der Probe nicht gemessen 
werden, sondern bezieht sich auf Literatur- 
angaben. 

Die Messung des Wandwiderstandes erfolgte 
bei Schallintensitaten, bei denen der Strémungs- 
widerstand noch unabhangig von der Schall- 
schnelle war, also im linearen Bereich. Der 
Betrag des Wandwiderstandes bei der Resonanz- 
frequenz wurde als akustisch bestimmter Stré- 
mungswiderstand in Tabelle III eingetragen. Als 
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statisch bestimmter Strémungswiderstand der 
Proben wurde der mit der Strémungsapparatur 
gemessene und auf die Geschwindigkeit Null 
extrapolierte Strémungswiderstand in Tabelle III 
aufgefiihrt. In Abb.6 sind die gemessenen 
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zur Kanalwandung mit wachsender Frequenz 
immer steiler wird [7], [10], [11]. Es leuchtet ein, 
da®B eine ahnliche Erscheinung bei Kanalen mit 
stark wechselnden Querschnitten auftritt, wie es 
bei den porésen Steinen der Fall ist. 


Tabelle III 


Heraklith- 


Pie Tuffstein 


Schichtdicke 


Wandabstand 5 
Resonanzfrequenz 280 
£ d, statisch 35 
£ d, akustisch 62 


Wandwiderstande fiir zwei der untersuchten 
Baustoffproben aufgezeichnet. Der zugehérige 
statische, auf die Geschwindigkeit Null extra- 
polierte Strémungswiderstand der Proben ist 
durch Parallele zur Ordinatenachse eingetragen. 


statisch bestimmter Stromungswiderstand 


=2 


0 0,5 1,0 


Abb. 6. Wandwiderstandsdiagramm fiir zwei 1,5 em starke 
porése Bausteinproben. X ist der Imaginirteil, 
R der Realteil des Wandwiderstandes, Z der 
Kennwiderstand der Luft; 
@--- @ Schwerbeton in 9 em Wandabstand, 

~ Schlackenstein in 4 em Wandabstand. 


x 


Man erkennt aus Tabelle III und aus Abb. 6 
eine deutliche Diskrepanz zwischen dem statisch 
und akustisch bestimmten Strémungswiderstand, 
wobei der akustisch gemessene stets der héhere, 
maximal der doppelte Wert war. Es ist bereits 
von dem Rayleighschen Kanalmodell bekannt, 
dafS der akustisch wirksame Strémungswider- 
stand bei héheren Frequenzen hiéher wird als der 
statische, weil der Abfall der Geschwindigkeit 


Schlacken- 
stein 


Schlacken- 


Schwerbeton Risshedsus 


8 cm 

450 Hz 
21 gem? 5-1 
38 gem-2s-1 


y) Messung des Strémungswiderstandes 

im nichtlinearen Bereich 

Die Abhangigkeit des statischen Strémungs- 
widerstandes der Bausteine von der Geschwindig- 
keit 146Bt vermuten, da®B bei diesen Steinen auch 
akustisch der Strémungswiderstand von der 
Geschwindigkeit, d. h. der Schallschnelle, abhan- 
gen wird. In diesem Falle wiirden sowohl der 
Wandwiderstand als auch der Schluckgrad der 
porésen Bausteine von der Schallintensitaét ab- 
hangen und daher akustische Nichtlinearitaten 
auftreten. In der Tat haben die folgenden Mes- 
sungen gezeigt, daB bei den untersuchten poré- 
sen Bausteinen eine Abhangigkeit des akustisch 
bestimmten Strémungswiderstandes von der 
Amplitude der Schallschnelle vorhanden ist. 

Die fiinf unter 8) gemessenen porésen Baustoff- 
proben wurden auch im nichtlinearen Gebiet 
untersucht. Akustisch geschah dies bei dem 
gleichen Wandabstand im Rohr wie im linearen 
Fall. Als Frequenz wurde méglichst genau die 
Resonanzfrequenz benutzt, d. h. der Wandwider- 
stand war annahernd eine rein reelle GréBe und 
daher gleich dem akustisch bestimmten Stré- 
mungswiderstand. Aus dem Pegelwert im Maxi- 
mum und Minimum sowie aus dem reellen Wand- 
widerstand konnte durch Berechnung die Ampli- 
tude der Schallschnelle unmittelbar vor der 
Probe gewonnen werden. Diese Schallschnelle 
entsprach der auBeren DurchfluBgeschwindig- 
keit bei der statischen Strémungswiderstands- 
messung. In den Abb. 7, 8 und 9 wurde der 
akustisch und statisch bestimmte Strémungs- 
widerstand in seiner Abhangigkeit von der 
Schnelle bzw. der DurchfluBgeschwindigkeit 
festgehalten. Der Effektivwert der Schnelle 
bzw. die Strémungsgeschwindigkeit ist als Ab- 
szisse aufgetragen. 

Man erkennt aus diesen Messungen, dal} die 
Nichtlinearitat des Strémungswiderstandes aku- 
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stisch im allgemeinen erst bei gréBeren Ge- 
schwindigkeiten einsetzte als statisch. Dies wird 
besonders deutlich in Abb. 7 bei der Probe aus 
Schlackenkiesbeton, bei der im akustischen Fall 
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fiihrte zu einem Anstieg des Strémungswiderstan- 
des mit der Geschwindigkeit, der ebenfalls etwa 
parallel zu dem Anstieg des statisch gemessenen 
Strémungswiderstandes verlief. 


100 110 120 125 dB 100 110 120 100 110 120dB 
100 T | 100 T T T T T : a 
: a 
ieee se a gems gems" akustisch 
80 7 == 
x 300 
akustisch h/ 
60 
Sfo- statisch 
‘ 200 
40 ee A r 
es we 
on 
a statisch te YA akustisch = | 
SE ard TTT 100 
20 20 =0}=-° 
og ane De 


0 b) ercniss 5 0 
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Abb. 7. Nichtlinearitit desStrémungs- 

widerstandes bei zwei porésen Bau- 
stoffproben der Dicke d; 
@ --- @ Heraklith, d = 2,5 cm; Mef- 
frequenz 280 Hz, x—x Schlacken- 
kiesbeton, d= 1,5 em; MeBfrequenz 
450 Hz. 


steinprobe; 


bis herauf zu 15 cm/s noch keine nennenswerte 
Nichtlinearitat auftrat, wahrend statisch sich 
schon lange eine, wenn auch nicht sehr starke 
Nichtlinearitat bemerkbar machte. Bei der 
Heraklithprobe in Abb. 7 war die Nichtlinearitat 
im statischen Falle starker als im akustischen. 
Bei einer Strémungsgeschwindigkeit bzw. Schall- 
schnelle von 3 cm/s waren der statische und aku- 
stische Strémungswiderstand gleich. 

Bei den anderen Baustoffproben kam im aku- 
stischen Fall nach einem deutlichen Ubergangs- 
gebiet bei héheren Geschwindigkeiten der Be- 
reich, in dem der Anstieg des Strémungswider- 
standes mit der Geschwindigkeit abnlich wie beim 
statischen Strémungswiderstand erfolgte, d. h. 
die Kurven fiir den akustisch und statisch ge- 
messenen Strémungswiderstand verliefen in ihrer 
Abhingigkeit von der Geschwindigkeit etwa 
parallel. Das akustische Ubergangsgebiet lag bei 
der Schlackenbetonprobe in Abb. 8 zwischen 
lund7 em/s, beidem Tuffbeton in Abb. 9 zwischen 
0,2 und 1 cm/s. Bei der Schwerbetonprobe in Abb.8 
war der Strémungswiderstand akustisch bis zu 
Schallschnellen von 0,5 cm/s unabhangig von der 
Schallschnelle. Bei héheren Schnellen setzte dann 
entsprechend Abb. 8 die Nichtlinearitat ein und 


Abb. 8. Nichtlinearitat desStrémungs- 
widerstandes bei einer 1,5 cm star- 
ken Schwerbeton- und Schlacken- 


Abb. 9. Nichtlinearitat desStrémungs- 
widerstandes bei einer 1,5 cm star- 
ken Tuffsteinprobe; Meffrequenz 
300 Hz. 


x—x Schlackenstein; MeBfrequenz 
450 Hz, @ --- @ Schwerbeton; MeB- 
frequenz 400 Hz. 


Die vorliegenden Untersuchungen wurden mit 
Unterstiitzung des Bundeswohnungsministeriums 
durchgefiihrt. Herrn Professor Dr. L. CREMER bin 
ich fiir die Anregung hierzu und wertvolle Rat- 
schlage sehr zu Dank verpflichtet. Ferner habe 
ich Herrn Professor Dr. H. Marer-LeErsnitz zu 
danken, daf ich die strémungstechnischen Unter- 
suchungen im Institut fiir Technische Physik der 
Technischen Hochschule Miinchen mit den dorti- 
gen Mitteln durchfiihren konnte. 


(Eingegangen am 25. Juli 1955.) 
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EXPERIMENTAL DETERMINATION OF THE DAMPING 
OF PULSATING AIR BUBBLES IN WATER 
IN THE FREQUENCY RANGE 100 to 300 ke/s 
by H. HaEsKke 
III. Physikalisches Institut der Universitat Gottingen 


Summary 


The bubbles investigated, having radii from 10 to 32 y, are placed in a water-filled, thick- 
walled glass tube with a diameter small in comparison with the wavelength in the water. The 
reflection coefficient of the bubble has been measured using a pulse method and from this the 
damping constant has been calculated. The investigations showed that the diameter and proper- 
ties of the measuring tube may influence the properties of the bubble. The tube walls cause images 
of the bubble outside the tube and the radiation resistance of the bubble is thereby altered in 
the same way as in a series of radiators. Therefore a correction factor for the radiation resistance 
has been calculated which allows a determination of the appropriate damping constant of a bubble. 
The results show that the damping of air bubbles in the frequency range 100 to 300 ke/s can be 
explained completely from the losses due to radiation and heat conduction. 


Sommaire 


Les bulles étudiées, ayant des rayons de 10 4 32 yp, ont été injectées dans un tube de verre a 
paroi épaisse, rempli d’eau, et dont le diamétre était petit devant la longueur d’onde dans l’eau. 
On a mesuré par une méthode d’impulsion le coefficient de réflexion des bulles, d’ot on a déduit, 
par un calcul, la constante d’amortissement. Les résultats ont montré que le diamétre et les carac- 
téristiques du tube de mesure peuvent avoir de influence sur les propriétés des bulles. La paroi 
du tube produit des images d’une bulle 4 l’extérieur du tube, de sorte que la résistance de rayonne- 
ment de la bulle est modifiée de la méme fagon que dans une série d’éléments rayonnants. On a 
done calculé un facteur de correction 4 appliquer @ la résistance de rayonnement, ce qui a permis 
de déterminer le coefficient d’amortissement approprié d’une bulle. Les résultats montrent qu’on 
peut expliquer parfaitement l’amortissement des bulles d’air dans la gamme 100 4 300 kHz au 
moyen des pertes par rayonnement et conductibilité calorifique. 


Zusammenfassung 


Die untersuchten Luftblasen mit Radien von 10 bis 32 u befanden sich in einem mit Wasser 
gefiillten, dickwandigen Glasrohr, dessen Durchmesser klein gegen die verwendete Wellenlange 
im Wasser war. Gemessen wurde der Reflexionsfaktor der Luftblase mit einem Impuls-Verfahren 
und daraus ihre Dampfungskonstante berechnet. Es zeigte sich, daB GréBe und Beschaffenheit 
des MeBrohres die Eigenschaften der Luftblase verandern kénnen: Durch Reflexion an den Rohr- 
wanden entstehen Spiegelblasen auBerhalb des Rohres; infolgedessen verandert sich der Strah- 
lungswiderstand der Luftblase, wie es beim Zusammenfiigen mehrerer Strahler zu einer Strahler- 
gruppe geschieht. Es wurde daher ein Korrekturfaktor fiir den Strahlungswiderstand berechnet; 
damit ist es méglich, die Dampfungskonstante einer Luftblase zu gewinnen. Aus den erhaltenen 
Werten geht hervor, da sich die Dampfung einer Luftblase im Frequenzbereich 100 bis 300 kHz 
vollstandig durch Abstrahlungs- und Warmeleitungs-Verluste erklaren 1aBt. 


1. Introduction 


An air bubble in water represents a system 
capable of oscillation, the relevant compliance 
being given by the compressibility of the enclosed 
gas and the mass being that of the surrounding 
water taking part in the oscillations. This com- 
pliance mass-system has a definite resonant fre- 
quency, which, according to MINNAERT [1], is 
given by the relation 


1 3x P 
fo= za V 0 (1) 


4 


as a function of the radius of the bubble ( fp = re- 


sonant frequency, Ro = radius of the bubble, 
K = C€p/cy of the enclosed gas, P = hydrostatic 
pressure, ¢ = density of the water). 

If exposed to a sound field, the bubble is 
excited to pulsations, the amplitude of which is 
determined exclusively by the damping if the 
bubble is in resonance. The damping is essentially 
caused by radiation and by thermal losses, i. e. 
a part of the energy withdrawn from the sound 
field is re-radiated in the form of a spherical wave 
and another part is transformed into heat during 
the polytropic compressions and expansions of 
the enclosed gas. These losses due to heat con- 
duction had been calculated by Prriem [2]. 
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In earlier papers by Meyer and Tamm [3], 
Laver [4], Exner [5] and Exner and Hampe [6] 
the damping constants in the frequency range 1 to 
100 ke/s were measured. The aim of the present 
work was the determination of the damping con- 
stants of air bubbles with still higher resonant 
frequencies. The measurements have been extend- 
ed to frequencies as high as the increasing experi- 
mental difficulties in handling single bubbles 
and the increasing inaccuracy allowed (up to 
300 ke/s). 

At the same time the underlying physical 
principles of the methods used here and in earlier 
investigations [5], [6] have been reconsidered as 
regards the influence of the measuring system on 
the bubble and its properties. For this, the 
measurements were made on bubbles in measuring 
tubes having different diameter and virtually 
rigid walls. 


2. Description of the measuring method 


The radii of the air bubbles investigated, 
having resonant frequencies between 100 and 300 
ke/s, are between 32 and 10u. Due to their small 
volume their rate of rise is very small; because 
of the increasing ratio of surface to volume the 
bubbles dissolve more quickly in the surrounding 
water the higher their resonant frequency!. Thus 
the lifetime of a 300 kc/s bubble is only about 
1 second. 

In view of these properties of the bubbles a 
pulse reflection method similar to that used by 
Exner and Hamre and adapted to the demands 
of the present work proved most suitable. 


a) Apparatus 


Fig. 1 shows the whole apparatus schemati- 
cally. The pulse generator I generates pulsed 
oscillations with a repetition rate of about 300 s~! 
which are fed to the transmitter S which is 
coupled to the measuring tube R through its 
pressure chamber. At some distance from the 
transmitter the probe of a microphone projects 
into the measuring tube, and the microphone 
output is recorded on the screen of an oscillo- 
scope. The sweep of the oscilloscope is synchro- 
nized with the repetition frequency of the pulses 
so that a stationary picture is observed as shown 
in Fig.2. The pulse radiated from the trans- 
mitter proceeds along the tube, generates the 
deflection (I) on the screen when passing the 
microphone, travels on, and is partly reflected at 


1 The water must be slightly undersaturated with air to 
avoid the occurrence of disturbing additional air bubbles. 
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Fig. 1. Schematic diagram of the 
apparatus (R measuring 
tube, S transmitter, M 
microphone, B air bubble, 
I pulse generator, A relay 
attenuator). 


the branch of the measuring tube and at its end, 
yielding another deflection (III). Further (second) 
deflections at the transmitter are not observed 
on the screen since at 
these high frequencies 
the damping in the 
measuring tubes is com- 
paratively high and is 
moreover partly in- 
creased by wax on the 
outside of the tube. 

If now the air bubble 
to be investigated (B in 
Fig. 1) is placed in the 
measuring tube the 
sound wave re-radiated 
from the bubble gener- 
ates an additional de- 
flection (II) on the 
screen. 


Fig. 2. Diagram of the pulses 
as recorded on the screen 
of the oscilloscope. 


b) Measuring procedure 


A water jet is directed momentarily into the 
water surface (at a in Fig. 3), thus generating a 
number of air bubbles of different diameter. 
These bubbles pass through the side branch of 
the cylinder b in which they circle and rise, so 
that it is easily possible to move some bubbles of 
the desired diameter into the perpendicular tube 
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c by a short pressure pulse. In this tube the 
bubbles rise more or less quickly until they collect 
in the horizontal part of the tube d arranged in 
order of diameter like a string of pearls. From 
among these bubbles a single one is chosen and 
moved by a current of water to the measuring 
point B. There it is kept motionless by means 
of a controllable stream of water directed down- 
wards, compensating its buoyancy. 

At first the selected bubble is somewhat too 
large to be excited to resonant oscillations by the 
chosen measuring frequency. During the decrease 
of the diameter of the bubble as it dissolves the 


Sapo 


— 


ne 
Fig. 3. Liquid-filled part of the apparatus made of glass 
(R measuring tube, a---f generation and selection 


of the bubbles, 4 bulb for generation of compressed 
air, 1---3 water circulation). 
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magnitude of the deflection (II in Fig. 2) on the 
screen of the oscilloscope; showing the pressure 
amplitude of the pulse reflected from the bubble, 
traces out a resonance curve as a function of 
time, the height of which is used for the caleula- 
tion of the damping constant. 

For this calculation, however, a knowledge of 
the pressure amplitude p,, with which the air 
bubble is excited, is necessary. To determine a 
reference magnitude a large air bubble several 
ems in length, filling the whole cross section of the 
tube is substituted for the bubble under test. 
Thus we get a deflection (II) which corresponds 
to total reflection at the position of the bubble?. 

For the calculation the diameter of the bubble 
in resonance must also be known. It is true that 
this diameter can be calculated from MINNAERT’s 
relation (eq. (1)) but to confirm this relation 
experimentally the diameter was measured by a 
second observer? at the same time as the measure- 
ment of the amplitude of the reflected pulse. For 
this purpose an image of the air bubble was pro- 
jected perpendicularly to the axis of the tube 
onto a screen with a 700-fold magnification, 


c) Technical details of the apparatus 


The pulse generator generated pulsed oscilla- 
tions of the measuring frequency fo. The pulse 
length was variable and was adjusted in such a 
way that the air bubble was excited to at least 
30 oscillations in the steady state. 

As transducer a barium titanate disk was 
used. It was coupled to the liquid within the 
measuring tube by means of a pressure chamber 
for which the exponential form proved most 
suitable to avoid resonances. 

The microphone used was a_ hydrophone 
designed by Tamm [7], with its conical probe pro- 
jecting into the interior of the measuring tube. 
The diameter of the tip of the probe was only 
0.4 mm and therefore always small in comparison 
with the wavelength in the water. 

At higher frequencies the pressure amplitude 
of the reflected pulse was often very small (the 
reflection coefficient of the pulse was sometimes 
below 3 %) so that the deflection (II) at the screen 
of the oscilloscope was too small for exact mea- 
surements. This difficulty could be overcome by 
using a relay with an attenuator (A in Fig. 1), 
which diminished the large primary pulse (from 
Ia to I in Fig. 2) and allowed a higher ampli- 


fication. 


2 The curvature of the reflecting meniscus is without 
influence on the measurement, as was established in pre- 
liminary experiments. 

3 Measurements made by E. MunpRy. 
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Fig. 4. The glass apparatus (R measuring tube, S transmitter. 
M microphone, B air bubble, a---f generation and 
selection of the bubbles, 4 generation of compressed 
air, 1---3 water circulation). 


The occurence of so-called ‘“‘anomalous”’ bub- 
bles during the investigations of ExNER and 
Hamper, which was obviously connected with the 
presence of dust particles or other impurities, led 
to the construction of a glass apparatus (Figs. 3 
and 4) which allowed extremely clean experimental 
conditions. This apparatus was built in the form 
of a completely closed system of tubes fused to 
each other, in which the water and even the air 
were circulatory. Where a flexible connection was 
essential (coupling of transmitter and hydrophone, 
pressure bulb), polyethylene, a plastic as inert as 
possible, was used. 

The measuring tube R was a thick-walled glass 
tube with a wall thickness 2 to 3 mm. The walls 
of this tube were virtually rigid in comparison 
with the water; measurements of the sound velo- 
city in the tube showed only a deviation of about 
1% from the sound yelocity in the free liquid. 

The measurements were made in different tubes 
with inner radii of 1.033; 1.209; 1.39 and 2.00 mm 
resp. In the frequency range in question these 
radii are in accordance with the two limiting 
conditions: To allow only plane waves to proceed 
within the tube the radius must be smaller than 
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the 0.6 of the wavelength in water; the cross sec- 
tion of the tube must exceed the absorption 
cross-section [5] of the air bubble. 

The buoyancy of the air bubbles was compen- 
sated by a corresponding downward stream velo- 


city applied during the period of the measure- 


ments. For this purpose the water was allowed 
to flow from the supply vessel 1 to the vessel 2 or, 
during longer measurements, the water was 
circulated by a glass propeller 3. 

The generation of the bubbles as well as all 
other operations within the glass apparatus were 
made exclusively by air pressure, produced by 
pressing the bulb 4 in a suitable way and led 
through the tube system to the desired place. 

The whole measuring apparatus was tested by 
measuring the reflection coefficient of small 
nickel cylinders placed in the measuring tube 
instead of air bubbles. A comparison of the mea- 
sured reflection coefficient with the easily cal- 
culated one provided the desired check on the 
method. 


3. Evaluation of the damping from the 
measurements 


As already mentioned two quantities were 
measured for the air bubbles investigated: 
the pressure amplitude p, of the exciting wave 
and 
the pressure amplitude p, of the re-radiating 
wave directed towards the transmitter. 
From these two quantities the damping con- 
stant4 of the air bubble can be calculated. 
EXNER [5] and Exner and HAmpE [6] used the 
following expression for the damping constant® 


Peet hai le py 
01 — Op + 030 — J) ( = fe (2) 
: e€ 


(6, and dso = damping constant due to thermal 
and radiation losses resp., 9c — characteristic 
impedance of the water, Rg = resonant radius of 
the bubble, r = radius of the measuring tube, 
fo = measuring frequency). 

For the calculation of the necessary value of 
dso they used the expression for the radiation 
resistance of a zero-order radiator with a dia- 
meter small in comparison with the wavelength: 
Zs0 = 4 Ri oc k, Ri and obtained 559 = Zs0/aM 
= 2zfo Ro/e. 


4 The damping constant is defined by “logarithmic 
decrement devided by x”. 

5 The values for the damping constants obtained in 
different ways are distinguished by the indices 1, 2 or 3, 
while the values finally resulting from the present investi- 
gation have the index 0 (see page (274)). 
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Fig. 5. Damping constants 6; as a function of the resonant frequencies of the bubbles derived by the method of EXNER 
and Hampr. Full-line curve calculated from radiation and heat conduction. The radius of the measuring tube 
was (a) 1.033 mm, (b) 1.209 mm, (c) 1.39 mm and (d) 2.00 mm. 


Applying this method to the new results we 
would obtain the values of the damping constant 
marked in Fig. 5a---d by points: 

These values depend on the radius of the meas- 
uring tube used and are moreover essentially 
smaller than the theoretical values in a certain 
frequency range. 

For comparison it was desirable to determine 
the damping constants also in another, direct 
way. For this purpose the sound pressure of the 
reflected pulses was measured with some bubbles 
as a function of the radius R of the bubbles 


(Fig. 6a). In Fig. 6b, this curve has been corrected 
by the factor (Ro/R)? and the damping constant 
62 can be taken directly from the relative half- 
width of the resulting resonance curve. This 
direct method yielded essentially higher damping 
constants 62 which, as is physically plausible, are 
even higher than the theoretical values. 
Therefore it seemed necessary to correct the 
relations used by EXNER and HAmMpPE to get a 
better description of the behaviour of the air 
bubble when in the measuring tube. For this 
purpose the calculation given above (eq. (2)) has 


30 A 30 


and it was found 
that the expression 
used for the radia- 
tion resistance need- 


been reconsidered 
] i | 
> 


25 


20 = — 


Pon 
> 


ed modification: 
The expression for 

the radiation resis- 

tance 

Zs0 = 4 7 R? ock2 R2 

is defined for a zero 


Fig. 6. Measured ratio of 
the reflected and 
incident 
pressures pr/p,e as 
a function of the 
radius R of the 


bubbles and same 
after correction 
by the factor 
(Ro/R)®. ~Meas- 
uring frequency 
fi .o= 136 ke/s, 
measuring tube c. 


sound — 
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order transmitter, situated in an infinite medium. 
During the measurements, however, the air bubble 
is in a narrow tube with rigid walls which cause 
images of the bubble outside the tube. Thus, 
when the bubble is brought into the measuring 
tube, its radiation resistance will be altered. 
For the calculation it is denoted Z; = Za, 
where a is a correction factor taking into account 
the influence of the tube walls. This factor will 
be calculated in a following section. 

The expression for the mass of the medium 
taking part in the oscillations M = 47 Ro, the 
value for the sound velocity ¢ in water and the 
relation of MINNAERT (eq. (1)), in a first approxi- 
mation, can be kept unaltered as calculations 
showed. Only the adiabatic exponent «x in eq. (1) 
will be replaced by the more exact polytropic 
exponent, in view of the high thermal losses in 
the frequency range under consideration. 


a) A new and more detailed description of the pro- 
perties of air bubbles 


The walls of the measuring tube are considered 
as rigid. Under the condition that the radius r is 
smaller than 0.6 wavelength, only a plane wave 
can proceed inside. The frequency of the incident 
wave may be denoted fo, the corresponding 
radius of the bubble in resonance Ro. Suppose 
a bubble with the arbitrary radius R has the 
resonant frequency f. 

If the bubble is excited into pulsating oscilla- 
tions by an incident plane wave with the pressure 
amplitude p, (the time-dependent factor ej 2*fot 
being always omitted) the surface velocity 
will be 


is Pe q 2. Peq (3) 
vA (Zy+ Zs) =} (wo M=1/a0 F) ’ 
where M is the oscillating mass, F the compli- 
ance, Z, the loss resistance, Z; the radiation 
resistance and q the surface area of the bubble. 
In unit time‘ the bubble radiates the energy: 


Up 


E, = = Zs|v» |? (4) 
ie ie |pel?q? Zs 
2 (Ly + Zs)? + (wo M—1/ F)? 
Since the walls of the tube were assumed to be 
rigid, the radiated energy can only proceed in 
both directions along the axis of the tube. The 
reflected wave directed towards the transmitter 
therefore has the energy: 


E, = E,/2 
and the energy of the incident plane wave (per 
unit time) can be expressed by 


& . em2 e s-3, 


6 The energy is here always referred to unit time. | 


H. HAESKE: DAMPING OF PULSATING AIR BUBBLES IN WATER 


_ 1 |pel? s 


2) (pe 


Ee (5) 


(S = mr? = cross-section of the tube). 
We thus have for the ratio between the ener- 


gies of the reflected and the incident waves, 


using the relation w? = 1/MF, 
LT fe 
z + ; 4 i, an ;2 (6) 
e (Zy + Zs)2 + w2 re -) 


wo 7) 


In the present experiments the frequency fo of 
the exciting wave was kept constant while the 
radius of the bubble and hence its resonant fre- 
quency f altered continuously as it dissolved. 
M, Z, and Z; therefore are not constants but 
functions of the radius R of the bubble: 


(i) M ~ 47 R80 as already given above. 


(ii) The loss resistance Z, includes all pro- 
cesses by which an irreversible transition of 
oscillatory energy into heat is caused. In the 
measurements with frequencies below 100 ke/s 
obviously only the losses by heat conduction as 
calculated by Prriem are of importance. It must 
be expected therefore that also in the frequency 
range 100 to 300 ke/s these thermal losses contrib- 
ute most of all. According to Prriem, in this 
frequency range 6; is nearly independent of the 
radius in the bubble. If we therefore assume 6;), 
to be independent of the radius—as a first 
approximation the same is valid then for dy, too— 
we get as an approximation for the loss resistance 
Z, the expression 

2 
Zy = by|oF oy ome (Ro fo)? R, 
where F=1/12 2 0g Cy R is the compliance, 0, 
the density, c, the sound velocity in air. 


(iii) The radiation resistance is expressed by 


Zs = Zso x, where the factor « is, to a first 
approximation, independent of R, as shown in 
the next section. 

Introducing these expressions for M, Z,, Zs, 
and using the relation Rofp = Rf = const (eq. 
(1)) and the definition of the damping constant 
6 = Z/wM we get 


mens) (7) 
x 725 x (R/Ro)* Serer 
[5 (Ro/R) + ds0 x (R/Ro)2]2 + my 


Graphic representation of this expression is 
impossible because of the unknown damping 
constants 6, and ds9. A good approximation, how- 
ever, can be achieved with the two assumptions: 
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(i) In the vicinity of the resonance, i.e. for 
R/Ro ~ 1 we have 


Oy = Os0 4 (R/ Ro)? i) dy ate Oso x 


for ds9 « smaller than dy. 
(ii) In the regions remote from resonance, 1. e. 
for R/Ro > 1 or <1 the inequality 
[dy + dso x (R/Ro)?]? <« [(R/ Ro)? — 1) 
holds. If we moreover introduce p,/p, = VE,|Ee 
and for brevity 63 = 6, + ds0 x we finally get for 


the dependence of the sound pressure of the reflect- 
ed wave on the radius R the desired expression 


Pr ~ V2x Ro (R/Ro)? 


f - Rad ES, 
Pe r Voi+ (R/Ro?—1P 


This quantity is plotted in Fig. 7 as a function 
of Ro/R with different values of 63 as parameter. 
The experimental points are taken from Fig. 6. 
For this, the factor « has been determined 
empirically by comparison of (p,;/P¢)caiculated 
with (p,/Pe)measured - 


Fig. 7. Ratio of the sound pressures p;/pe as a function of 
the radius R of the bubbles calculated according 
eq. (8). The experimental points are taken from Fig. 
6a, measuring frequency fo = 136 ke/s. 


Fig. 8 shows the same quantity for another 
frequency. 

The good agreement of the calculated curves 
with the experimental results allows the conclu- 
sion that the calculation suitably describes the 
behaviour of the bubble. If further the resonance 
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Fig. 8. Experimental results and calculated curve for p;/pe 
as a function of the radius of the bubbles for the 
measuring frequency fo = 221 ke/s. 


maximum of the experimental curve in Fig. 7 is 
compared with the calculated curves, this com- 
parison yields the same value (6 = 0.19) of the 
damping constant 63 as the one formerly taken 
directly from the measurements in the form of 
a relative half-width (Fig. 6)—a further con- 
firmation of the validity of the above considera- 
tions. 

Summarizing, we can describe the influence of 
the tube wall on the physical properties of the 
air bubbles as follows: 

The radiation resistance of the bubble to be 
investigated is altered by the existence of image- 
bubbles in the same way as in a series of radia- 
tors. The energy radiated by the bubble is 
altered correspondingly, the sound pressure being 
larger (or smaller) by the factor )/a- 

As a consequence of this altered radiation the 
air bubble in the tube really has a damping 
constant different from that effective in the 
infinite medium, namely 62 = 43 = 6, + ds0 x. 


b) Calculation of the correction factor for the 


radiation resistance 


According to the above considerations a correct 
description of the behaviour of an air bubble in 
the measuring tube is possible by means of a 
correction factor x. Attempts have therefore been 
made to calculate this factor: 
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x) One dimensional case 

The image-bubble model is not helpful for the 
circular measuring tube. As a first step therefore 
the one-dimensional simplification was considered. 
The air bubble was assumed to be situated be- 
tween two infinitely extended rigid planes, at a 
distance 2r from each other. It will be assumed 
that r is large in comparison with the radius R 
of the bubble. The equivalent model, an infinitely 
extended line of discrete radiators of order zero, 
placed in an infinitely extended medium, has 
been taken as a basis for the following calcula- 
tions: 

All bubbles of the line are excited to cophasal 
oscillations by a plane sound wave incident 
normally to the line. Each of them then radiates 
a spherical wave with the sound pressure 


P(r) = (A/r) 3" 
and the sound velocity 
v(r) = (A/oer) (1 —j/kor) e*. 
At the position of one bubble the contributions 


to the sound pressure from all radiators add 
up to 


Pw = P(R) + 2S. p2rp)e%". (9) 


p(R) is the sound pressure generated by the 
bubble under consideration; the factor e~* takes 
account of losses due to the incomplete reflection 
at the tube walls. 

The energy withdrawn from the bubble per 


unit time by this total sound pressure equals: 
N= 5 qe(R)* Pro 5 (10) 


the acoustical impedance Z, of a radiator in 
general is defined by the relation: 


N= + Z, |o(R)[2, 


where N is the sum of the radiated energy and 
the reactive power. 

By equating these two expressions (eqs. (10) 
and (11)) one gets for kgR < 1 the impedance 


Zs = qock, R2(1+ j/ko R) 


f ae >> i ositn,| : (12) 
rom 
= 


Caleulation of the sum for 0 < kor < x yields 


(11) 


ps =—jln (1 + e-26—e-8 cos 2 kor) — 
e-8 sin 2 kor 
1—e-*cos2 kor ~ 


—j arctan - 


In the limiting case of complete reflection at the 
tube walls, i.e. for 6 +0, one gets instead the 
_ expression 
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™ — 2kor 


lim [5°] = —In[2 sin kor| —j 


p—-0 

(the evaluation of which is easier) and hence we 
get finally for the radiation resistance of an air 
bubble between two reflecting planes 


Zs= Re(Z;) 


Ais daly coe. _n—2kor 
= qgek; R (1——"In)2 sin kor | + j ana 
= ZL 50 x] (0 <kor 7). (13) 


On account of the assumption introduced that 
r > R which was always fulfilled in the experi- 
ments we can neglect the second term of the 
correction factor. Thus «; is independent of the 
radius of the bubble, as was assumed above. 
1 —2kor 
2 kor 


In Fig. 9 this correction factor «; ~ 1 + 


is represented as a function of kor. 


si ecaars 


Fig. 9. Correction factors for the radiation resistance as 
functions of frequency and radius of the tube. «,; has 
been calculated for one-dimensional images, «2; and 
&s are approximate solutions for the conditions in 
the measuring tube. The values marked by circles 
have been determined from measurements. 


6) Two-dimensional case 

The calculation of the correction factor for the 
case of images in two dimensions was difficult. 
An attempt was made to treat the problem of 
an air bubble in a circular tube as well as in a 
square one (both problems ought to be physi- 
cally similar) but in both cases the calculations 
led to sums which cannot be expressed in a 
closed form, so that only a numerical evaluation 
for some special cases would be possible. 

If, however, the procedure is applied of multi- 
plying the one-dimensional direction factors with 
each other to get the directivity of the two- 
dimensional arrangement, which is exactly per- 
missible in the case of diffraction calculations 
referring to some distant point, the above cal- 
culated one-dimensional correction factor «; leads 
to an approximately valid solution for the pro- 
blem of the air bubble in the measuring tube. 


274 


This factor «; is also represented in Fig. 9 as 
a function of kor. 

A further approximate solution for the round 
tube can be got in the following way: With a 
tube of square cross-section the images would be 
distributed over the whole plane. With a round 
tube the images form concentric rings having the 
constant distance 2r from each other. A descrip- 
tion of the behaviour of these ring-shaped bubble 
formations is difficult, but since there cannot be 
any serious difference in the physical effects of 
both the arrangements, a constant number of 
bubbles per unit area over the whole plane may 
be formally introduced for the case of the round 
tube, too. The sound pressure contributions of 
these image bubbles are summed now over 
circular rings of the width 2r and we get the 
expression 

&=1+ Dy ree 

kor 
This rough approximation to the correction fac- 
tor has also been represented in Fig.9 for 
comparison. 

In the same figure moreover some empirical 
values for « are given which have been deter- 
mined by comparison between the measured and 
calculated minima of the curve for p;/pe. These 
values confirm the usefulness of the approxima- 
tion a;; their deviations from the approximation 
curve are systematic and are caused by the width 
of the lines on the oscilloscope screen which was 
especially disturbing in measurements of the 
minima. 


(0 <kor <7/2). (14) 


4. Results of the damping measurements 


With a suitable description of the behaviour of 
the bubbles in the measuring tube and with the 
necessary correction factor « approximately 
known, the desired damping constants could then 
be calculated from the measuring results. 

Corresponding to eq. (2) one gets from eq. (8) 
for the damping constant of an air bubble in an 
infinitely extended medium 


: —R 
do cs Oy+ 030 = V2 (o4 : Be Os0 (a A): 


(15) 
r pr 

The mean values of the damping constants 
given in Fig. 5 a---d which had been determined 
with the method of ExNER and Hampk (eq. (2)) 
have been corrected according to eq. (15) using 
the approximate value «7. The damping con- 
stants for an air bubble in an infinite medium 
got in this way are represented in fig. 10: 

The influence of the radius of the measuring 
tube has vanished, the values resulting from 
measurements in the different tubes a to d agree 
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Fig. 10. Damping constants 69 of air bubbles in water as 
a function of their resonant frequencies. (The 
curves a---d have been derived from mean values 
of the results given in fig. 5 a---d by correction with 


a? The full-line curve gives the sum of the cal- 


culated damping due to radiation and heat con- 
duction.) 


well with each other and prove the suitability of 
the correction factor. Further confirmation is 
given by those experimental points for which the 
correction factor is unity and which therefore 
could be taken directly from Fig. 5 a---d. These 
values have been marked by squares in Fig. 10. 

The experimental values for the damping 
constant are only slightly higher than the full 
line giving the sum of damping by radiation and 
damping by heat conduction calculated according 
to Prriem. Obviously therefore the damping of 
a pulsating air bubble in water is also determined 
in the frequency range 100 to 300 ke/s by losses 
due to radiation and heat conduction. Other 
causes of damping possibly effective, as for 
example the viscosity of the surrounding water 
or a possible hysteresis of the surface tension, do 
not count in comparison with the large contri- 
bution of the thermal losses. Their contributions 
are at least smaller than the possible error due 
to the inaccuracy of measurement. 


5. Retrospective consideration of earlier investi- 
gations on air bubbles in water 


At the end of this paper something may be 
said concerning the measuring method used in 
earlier investigations, especially these of EXNER 
and Hamper [5], [6] in which the necessity of a 
correction factor did not appear: 

The correction factor introduced during the 
present investigations is valid for totally refleet- 
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ing walls of the tube, i.e. for a tube which is 
absolutely rigid. The other limit, a tube in which 
no images occur at all, would be given by a tube 
the walls of which have an input impedance equal 
to the characteristic impedance of the water. 
The measuring tubes used by the various authors 
had properties between these two limits. 

Thus: 

In the present investigations the tubes used 
were thick-walled and little excited to vibration, 
nearly realizing the case of totally reflecting walls; 
EXNER and Hampsr, however, used tubes with 
a smaller ratio of wall thickness to tube diameter. 
Thus, their tubes took part in oscillation to a 
higher degree, resulting in a smaller influence of 
the images. 

(The decrease of the sound velocity in the tubes 
used here against that of the free medium is only 
1% against 4% in EXNER’s and even about 10% 
in Hampe’s tubes.) 

A further remark may be added: as already 
described the system of glass tubes used for the 
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present measurements was completely closed, so 
that the purity of the liquid could be kept much 
higher than in the former experiments. The 
occurrence of so-called ‘‘anomalous’”’ bubbles 
(i.e. bubbles with resonant frequencies higher 


than those expected from eq. (1) [6]), has, per- 


haps as a consequence of that, never been 
observed. 


The author wishes to thank Prof. Dr. E. MEYER 
for his kind interest and for many suggestions, 
and is further indebted to the Department of 
Scientific and Industrial Research, London, for 
having supported this work. 
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KOLLOQUIUM FUR AKUSTIK UND ULTRASCHALL IN UNGARN 


Die Physikalische Gesellschaft ,,L6rand Eétvés“ rief 
vom 21. bis 23. November 1955 in Budapest ein Kolloquium 
fiir Akustik und Ultraschall zusammen, auf dem sich die 
Fachleute dieses Gebietes zum ersten Mal zum Meinungs- 
austausch trafen. Die sieben zusammenfassenden V ortrage 
und 20 Berichte iitber Forschungsergebnisse, die zum Vor- 
trag kamen, wurden rege und eingehend diskutiert, auch 
von Seiten der ZuhGrer. Sowohl die etwa 33 Teilnehmer der 
Tagung, als auch die Vertreter wissenschaftlicher Vereine 
und der Ministerien waren~sich iiber den groBen Erfolg 
einig und gaben der Meinung Ausdruck, dafs damit fiir die 
Zukunft der Grund fiir eine koordinierte Arbeit gelegt sei. 
Es wurde beschlossen, alljahrlich ahnliche Referaten- 
tagungen abzuhalten. 

Auf dem Gebiet der Akustik wurde iiber Fragen der 
Psychoakustik, der Raumakustik und der Lautsprecher- 
entwicklung berichtet, auf dem Ultraschallgebiet iiber 
Materialuntersuchungsprobleme und zahlreiche Ergebnisse 
der physikalischen und physikochemischen Ultraschall- 
wirkungen. 

Die Themen der Vortraige waren: 

Tarnoczy, T., Die neuesten internationalen Ergebnisse der 
Sprachlautforschung (zusammenfassender Vortrag). 
Misxouczy-Fopor, F., Zur Psychophysik der Lautheit (zu- 

sammenfassender Vortrag). 

Macyart, E., Akustische Gesichtspunkte bei der Klassi- 
fizierung der Geraiusche (zusammenfassender Vortrag). 
Bresoyszky, J., Widerstandsinderungen des Mikrophon- 

Graphits bei Einwirkung von Gleich- und Wechseldruck. 
Vatk6, I. P., Anwendung elektroakustischer Gesichts- 

punkte bei der Untersuchung der Réhrenmikrophonie. 
Bardr, Z., Wirkungsgrad des direktstrahlenden dynami- 

schen Lautsprechers. 


Huszty, D., Annihernde Theorie des in einer unendlichen 
Wand strahlenden Ringes. 

Ferner, E., Erfahrungen und Messungen mit heimischen 
Druckkammerlautsprechern. 

Macyari, E., Akustische Regelungen in Hérsalen durch 
Resonatoren. 

SeB6K, L., Bedeutung der Geriuschpegelmessungen (zu- 
sammenfassender Vortrag). 

Luxics, M., Zur Erhéhung der Genauigkeit der Absorp- 
tionsmessungen nach dem Hallraumverfahren. 

ZspDANSZKY, K. und Fonacy, I., Ein Tonhdhenschreiber 
und einige mit diesem durchgefiihrte Untersuchungen. 

KeEMENES, F., Verstandlichkeitsuntersuchungen in Licht- 
spieltheatern. 

Ketter, G., Untersuchung der Schwingungsdimpfung von 
GuBeisen mit Ultraschall. 

Eert, [., MeBtechnische Probleme bei der Immersions- 
Durchleuchtung von Stahlplatten mit Ultraschall. 

Tarnoczy, T. und Invényt, A., Anderung der Schallge- 
schwindigkeit und des Schallwellenwiderstandes in einem 
Gemisch von Petroleum und Lein6él zwischen 10° und 50°C. 

SzmArp, J., Schallgeschwindigkeits- und Absorptions- 
untersuchungen an Beton mit dem Ultraschall-Echolot- 
gerat. 

Haspu, S., Messung der Durchgangszeit einmaliger Schall- 
impulse durch Beton. 

Retr, P., Méglichkeiten und Aussichten der Materialunter- 
suchungen mit Ultraschall (zusammenfassender Vortrag). 

Grecuss, P., Industrielle Bedeutung der Anwendung von 
Ultraschall in der Chemie (zusammenfassender Vortrag). 

Kamocsay, D., Physikalische Probleme bei der Anwendung 
des Ultraschalls in der Medizin (zusammenfassender 
Vortrag). 


Tarnoczy, T., Vom sogenannten Ultraschall-Kreuzfeld. 

Grecuss, P., Herstellung von Bitumen-Wasser-Emulsion 
durch Ultraschall. 

Rusznix, I., Grou, M. und Tarnoczy, T., Anderung des 
Dispersitatsgrades von Farbstoffen bei Einwirkung von 
Ultraschall. 

Tamas, Gy. und Tarnoczy, T., Verinderung der Durch- 
lassigkeit von Membranen bei Einwirkung von Ultra- 
schall. 

GouiAn, B. und Tamis, Gy., Auszug des Alkaloids aus 
Belladonna mit Ultraschall. 

DomAn, P., Verainderung des Redox-Potentials bei Ein- 
wirkung von Ultraschall. 

Als Gegenstand des Kolloquims 1956 wurden die Fragen 
der Schallempfindung, Geraiuschmessung und der Larm- 
bekampfung gewahlt. 

T. Tarnéczy, Budapest. 


Buchbesprechung 


Klangstruktur der Musik, Neue Erkenntnisse musik- 
elektronischer Forschung, im Auftrage des AuBen- 
instituts der Technischen Universitit Berlin-Charlotten- 
burg zusammengestellt und bearbeitet von F. WINCKEL. 
Verlag fiir Radio-Foto-Kinotechnik, Berlin-Borsigwalde, 
224 Seiten, DIN A4, 140 Abb., Ganzleinen DM 18,50. 
Das Buch enthalt eine Zusammenstellung der Vortrage, 

die im Jahre 1954 in der Vortragsreihe ,,Musik und Technik“ 

des AuS eninstituts der Technischen Universitat Berlin- 

Charlottenburg gehalten wurden. Durch drei weitere Bei- 

triage ist die Zusammenstellung ergaénzt worden. 

Man kann nicht erwarten, daf sich eine Vortragsreihe 

» Musik und Technik*‘ unter das Thema ,,Klangstruktur 

der Musik; Neue Erkenntnisse musik-elektronischer For- 

schung® subsumieren lat. Tatsaichlich wird in dem Buch 

z. B. nichts tiber die speziellen Strukturen der Klange der 

verschiedenen herkémmlichen Musikinstrumente mit Aus- 

nahme der Orgeln gesagt. Andererseits bringt das Buch — 
und nicht zu seinem Schaden — sehr viel mehr, als der 

Titel verspricht. Eine Angabe der Beitriige mége das demon- 

strieren: 

I Wincxke.L, F., Naturwissenschaftliche 
der Musik. 

Il DrAcer, H.-H., Die historische Entwicklung des 
Instrumentenbaues. 

LorrermoserR, W., Akustische Untersuchungen 

an alten und neuen Orgeln. 

IV Saua, O., Subharmonische elektrische Klangsyn- 

thesen. 
V Pouttin, J., Musique Concréte. 

VI Meyer-Eppier, W., Elektronische Musik. 

ENKEL, F., Die Technik des Tonstudios. 

STEINHAUSEN, H.-W., Musische Technik. 

IX Bracuer, B., Die musikalische Komposition unter 

dem EinfluB der technischen Entwicklung der 
Musik. 

X Sruckenscumipt, H.-H., Musik und Technik. 
(Leider ist auch aus diesen Angaben nicht immer klar zu 
erkennen, was im einzelnen behandelt wird. Das jedoch 
naher anzugeben, fiihrte hier zu weit.) 

Die Autoren sind Ingenieure, Musiker, Musikwissen- 
schaftler und Physiker. Entsprechend sind die Beitrage in 
ihrer Art recht verschieden. Herausgeber und Verlag haben 
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das Buch fiir Musikwissenschaftler, Komponisten, Musiker, 
Ingenieure, Tonmeister, Physiker, Elektroakustiker, Musik- 
instrumentenbauer und interessierte Laien bestimmt. Und 
gewiB werden alle einen Gewinn davon haben. 
Fiir die Leser der Acustica kénnen die Beitrage I, IIT, IV, 
V und Vials Zusammenfassungen (mit Literaturnachweisen), 
II, VIII, IX und X fiir die ,,akustische Allgemeinbildung~ 
wertyoll sein. Der Beitrag VII enthalt in seinem 2. und 
3. Teil leider einige Unrichtigkeiten und Unklarheiten, die 
in einer 2. Auflage des Buches vermieden werden sollten. 
G. R. Schodder. 


Heinrich Barkhausen T 


Am 20. Februar 1956 verstarb in Dresden Prof. Dr. phil. 
Dr.-Ing. E.h. Heinrich Barkhausen im Alter von 74 Jahren. 
Mit ihm verliert die Wissenschaft einen der groBen Forscher 
und Lehrer, die durch ihre eigenen Arbeiten und durch ihre 
Schule weit iiber ihr Fachgebiet hinaus zu Ruf und Ansehen 
gekommen sind. 


H. Barkhausen wurde am 2. Dezember 1881 in Bremen 
geboren, studierte Physik in Berlin, Miinchen und Géttin- 
gen und promovierte 1907 bei H. Simon in Géttingen mit 
einer Arbeit iiber schnelle elektrische Schwingungen, die 
groBes Aufsehen erregte. Er arbeitete dann einige Jahre in 
den Siemens-Laboratorien und iibernahm im Jahre 1911 
den Lehrstuhl fiir Schwachstromtechnik an der Technischen 
Hochschule Dresden, wo er das erste schwachstromtech- 
nische Institut Deutschlands griindete, und wo er bis zu 
seiner Emeritierung tatig war. 

Seine Arbeiten haben seinen Namen untrennbar mit der 
Entwicklung der Schwachstromtechnik vyerkniipft, zu 
deren Vatern er gehért, doch ist hier nicht der Ort, Bark- 
hausens Verdienste um die Schwachstromtechnik und ins- 
besondere um die Entwicklung der Elektronenréhren zu 
wiirdigen. Die Akustik dankt ihm unter anderem die 
Methode der subjektiven Lautstirkemessung und die Ein- 
fiihrung des ,,phon* als Einheit der Lautstiarke. 

Barkhausens umfassendes Lebenswerk fand iiberall 
Anerkennung, wie die groBe Zahl internationaler Ehrungen 
und Auszeichnungen zeigt, die ihm zuteil wurden. 
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